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Preface 


This publication contains the presentations made at the Industry-wide 
Workshop on Computational Turbulence Modeling, which was hosted by 
ICOMP/LeRC, and took place on October 6-7, 1994 at the Ohio Aerospace 
Institute. The purpose of the workshop was to initiate the transfer of technol- 
ogy developed at Lewis Research Center (LeRC) to industry and to discuss 
the current status and the future needs of turbulence models in industrial 
CED. To address the latter, a total of fourteen presentations were made by 
researchers from industry. CMOTT would like to thank all the workshop 
speakers for bringing to our attention a host of problems which are impor- 
tant to industry and for which they think CMOTT can be of help. We are 
prioritizing all the suggestions in order to incorporate them into the CMOTT 
work plan. 

One unanimous recommendation of the workshop participants was to 
make the workshop an annual event. This first workshop grew out of the rec- 
ommendations by the peer review committee of the LeRC turbulence mod- 
eling program, held in September of 1993. It could have not successfully 
transpired without the help and guidance of Dr. Chander Prakash (GE- 
Aircraft Engines), Dr. Munir Sindir (RocketDyne), and Dr. Saadat Syed 
(Pratt & Whitney), and for this CMOTT would like to thank them. 
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TURBULENCE PROGRAM FOR PROPULSION SYSTEMS 
Tsan-Hsing Shih 

Institute for Computational Mechanics in Propulsion and 
Center for Modeling of Turbulence and Transition 
NASA Lewis Research Center 
Cleveland, Ohio 

BACKGROUND 

• CMOTT group at LeRC has been in existence for about 4 years. In the 
first 3 years, its main activities were in developing and validating tur- 
bulence and combustion models for propulsion systems, in am effort to 
remove the deficiencies of the existing models. Two workshops on com- 
putational turbulence modeling were held at LeRC (1991, 1993). 

• A peer review of turbulence modeling activities at LeRC was held in 
September, 1993. Seven peers (GE, P&W, RocketDyne, Cornell, Berke- 
ley and NASA Ames) conducted the peer review. The objective of the 
peer review was to assess the turbulence program at LeRC/CMOTT 
and to suggest the future direction of turbulence modeling activities for 
propulsion systems. 

• Important messages from the peer review: 

0 “LeRC should spend substantial effort being responsive to indus- 
try’s current pressing perceived needs; this involves extensive dis- 
cussion with industry during every phase of model development, 
analysis of industry’s problems, goal oriented model development, 
evaluation of models relative to industry’s intended application 

0 “LeRC has an obligation not only to respond to industry’s re- 
quests for help, but to play an autonomous, independent leader- 
ship role in providing models of the highest quality, . .,. which can 
be employed not only by the aircraft gas turbine and rocket indus- 
tries but also by other industries 

<0* “In the present financial climate, industry does not have the re- 
sources to undertake model development and evaluation. LeRC’s 
help in this regard via the creation of its turbulence modeling ef- 
fort, is, therefore, welcome from the industry’s standpoint.” 

“It is important to work with the industry to evaluate the models 
and rank-order them by performance and cost in order to identify 
the most appropriate models for particular situations.” 

0 Many other useful suggestions and comments including collabora- 
tion with industry, joint programs, industry-wide workshop . . . 
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PROGRAM GOALS AT CMOTT 


• Develop reliable turbulence (including bypass transition) and combustion 
models for complex flows in propulsion systems 

• Integrate developed models into deliverable CFD tools for propulsion 
systems in collaboration with industry. 

PROGRAM APPROACH 

• Develop turbulence and combustion modules for industry customers 

• Industry collaboration and technology transfer 

• Model development for propulsion systems 

0 One-point moment closures for non-reacting flows 
<£> Scalar PDF method for turbulent reacting flows 
0 Vali dation of existing and newly developed models 
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Development of 

Turbulence and Combustion Modules 


• Objective 

0 Build a quick and efficient vehicle for technology transfer to indus- 
try 

• The features of the turbulence module: 

<> It contains various turbulence models from which users can choose 
the appropriate model for flows of interest 

0 It is self-contained, i.e., it contains its own solver for turbulence 
model equations 

<£ It can be easily linked to industry’s CFD codes 

• Turbulence module for NPARC code has been developed, tested, and is 
ready to be released 

^ The models built-in at the present time: 

Mixing length, Chien k — e, CMOTT k — e models 
<£ The model to be built-in: 

CMOTT algebraic Reynolds stress, Reynolds stress transport 
equation models and other models based on the request from in- 
dustries. 

o Built-in robust, realizable numerical solver for model equations. 

• General turbulence modules 

0 Can be used for both compressible and incompressible flows. 

0 Interface programs for different industry CFD codes 
<> Built-in models will be periodically updated. 
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Nozzle exit height = 0.122 in. 



Total Temperature at x=10.5 in 
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Collaboration with Industry 
and Technology Transfer 


• Joint research pro gr aims with industry 

0 Preliminary programs with engine companies and others have been 
initiated (GE, P&W, RocketDyne, Naval Research Laboratories) 

0 Develop further joint research programs related to the industry’s 
projects 

• Industry-wide workshops will be a regular program (once every two 
years) 

0 Release Lewis turbulence and combustion modules to industries 
0 Discuss the needs of industry 

Models developed at CMOTT 


1. Isotropic eddy viscosity models 

2. Reynolds stress & scalar flux algebraic equation models 

3. Second moment transport equation models 

4. Multiple-scale models for compressible turbulent flows 

5. Bypass transition models 

6. PDF models for turbulent reacting flows 


PROGRAM SUMMARY 



6 











Isotropic eddy viscosity models 


• Objective 

^ To examine the deficiencies of existing models 
To develop better eddy viscosity models 

• Current status of existing k — e eddy viscosity models 

2 k 2 

—UiUj = vx{Ui y j + Uj,i) — - k8ij , i 't — 

— = T™ + P (fc) -e + W.C., — = T* e) + P (e) - D (c) + W.C. 

Dt Dt 

<C> They are not tensorially invariant due to f^{y + ), W.C.(y + ) 

^ Model constants are not consistent for flows with and 
without wall 

<^> Normal stresses may violate realizability 

<> Do not work very well for flows with pressure gradients 

• Development of a Galilean-, tensorially invariant, realizable, k - e model 

§ New d am ping function fn{k/S v) is proposed to remove the depen- 
dence on y 

New dissipation e equation is introduced to give better 
response to pressure gradients 

Consistent model coefficients for all flows 

Realizability of the normal stresses is guaranteed 

Modified wall function for cases with pressure gradients 
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• CMOTT k — e eddy viscosity model 


2 k 2 

-muj = v T {Uij - Uj,i) - - k8ij , v T - Cpfp — 

O £ 


Dk 

Dt 


= Tk + Pk - £ 


De 

Dt 


— T e + CifiS e — C 2 





0 are functions of R = k/Sv , which is tensorially invariant 

^ ~ A^TaTU* fc/e’ which ensures realizability for normal stresses 

0 $ represents the effect of inhomogeneity 

k 2 k A 

# = Vk + b2 — VS Vk + 63 —VS VS 

£ £ 2 


• Validation 
Flows: 

0 Channel flows 

^ Boundary layer flows with and without pressure gradients 
^ Planar jet, round jet and mixing layer 
0 Backward-facing step flows 

^ Complex flows related to industrial applications 
Models: 

Q Launder- Sharma, Lam-Bremhorst, Chien, 

Nagano-Hishida, . . . 

$ k — u model (Wilcox) 

<> CMOTT k — £ model 
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Present model with the modified wall function 









Spreading Rate of Free Shear Flows 
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st. k — e 
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Algebraic Reynolds stress models 


• Objective 

^ To examine the deficiencies of existing ARS models 
^ To develop better ARS models 

• Current status of ARS models 

0 Second-order closure based ARS models (Rodi, 1980) 

—r-^-(P - e) = -TwTkUj'k ~ UjUkU it k ( p,iUj 4- p,jUi) - 2uv~Ju~j e 

k p 

Comments: 

* Assumption: UiUj/k = Const., (uiUjUk),k = (kui) t i = 0 

* Numerical difficulties 

^ Pope’s explicit ARS model (2-D flows), Taulbee’s ARS model (3- 
D), Gatski and Speziale’s ARS model 

0 Other methods: RNG, DIA and invariant theory 
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• General constitutive relations from invariant theory 

vZuj = \k8ij + 2a-2 — {Uij + Uj t i - - Ui,i8{j ) + 2a^—(U i j + U jti - - 
+ 2a 6 ^r(Ui >k U j ,k - ^n 2 6;j) + 2a 7 ^-{U k ,iU k ,j ~ -^ 2 hj) 


+ 2 a*^-{Ui, k Ul k + Uf >k Uj tk - |n 3 5ij) + Zaio—iUk'iUlj + U k<j Ul. 
+ 2 a 12 ^{Uf ik Ul k - \lU 6 ii) + ~ 1 ^) 

+ 2ai4^-(Ui, k Ui, k Uij 4- U i>k U ltk Ul - -Us**) 

Tf6 _ 2 V 

+ 2ai6-^(Ui t kUi tk Uij + U j>k U l<k Ui ti - -n 6 5ij) 

+ 2a ls ^-(U i<k U l , k Ul m Ul m + U jtk Ui, k Ul m Ul m - |n 7 6 0 ) 


• RDT and realizability constraints (Reynolds, Lumley) 

• CMOTT algebraic Reynolds stress model 

u jnj = I kSii - c„— 2S’. + 2 c 2 $(-s&nj y + n* t sy 

o £ ^ 


k t t + Ujkj = [(u H — ~)k,j],j - UiUjUij e. 

£ £ 2 
£,t + Uj£,j — \{v + — )e,j],j ~ C c \-UiUjUi,j - Cc 2 ~£ 


where 


C »- Ao + A;^' 


C 2 = 


y/l -QCli^y 

Co+ QSLhnik 


k? 


Vi 


= C„ — , Ao = 6 . 5 , Co = 1.0 


C c \ — 1.44, C e 2 — 1.92, <7fc — 1, o’c 1.3 


ThSij) 


i - §n»««) 
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• Validation 


0 Rotating homogeneous shear flows 
<C> Backward-facing step flows 
0 Confined jets 

0 Complex flows related to industrial applications 



Configuration of rotating Homogeneous shear flow 



Evolution of turbulent kinetic energy with time. 
— : present model; — - — : SKE; • : LES 
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Scalar turbulence model 

• Objective 

To improve the predictive capability of current scalar 
turbulence (0 2 - eg) models 

<£ A new scalar flux constitutive relation 

0 A new scalar dissipation rate model equation 



_ (2 + 2 r + 0.5r 2 ) 

A ~26 + 3.2t7 2 + 2£ 2 

St = v /0 .i 0 ,«> s = ^J2Sij Sij , Tj = Sk/e , £ = 7 (=) 1/2 Sr, = 

C« — Ci — 0.13, (7,92 = 0.63, C 93 = C 2 — 1, <7t = 1.0, = 1.8 
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Tavoularis & Corrsin Expt (1981) 



x/H 


Flat plate boundary layer with constant surface temperature 



y 
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• Validation 

0 Homogeneous turbulence subjected to dO/dy 

Q Homogeneous turbulence subjected to dU/dy, dQjdy 

0 Flat plate boundary layer with constant surface 
temperature 

• Work in progress 

§ Model assessment for different scalar boundary conditions 
^ Model extension for integration to the wall 


Sirivat & Warhaft Expt (1983) 
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Second Order Closure Models 


DuiUj 

Dt 


= T t j + p {j + n* apid + ng 


Return 



• Objective 

0 To assess existing models 

^ To find the direction of improving closure models 

• Basic model forms 

n**- d = 

Tij — F ij ((ujUj ),&, k, s') 

• General comments on second order closures: 

0 The model, Hg apid , is relatively well developed compared with 
other terms 

0 The model, ng eturn 5 is least developed 

^ A Galilean and tensorially invariant second order closure model 
has not been well developed yet 

0 All models have large errors near the wall, especially in the buffer 
layer; therefore, for engineering application, the wall function ap- 
proach is suggested at the present time 
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Application of realizability to IP and LRR models 





Multiple scale k-e model 


• Objective: 

0 To consider the effect of a non-equilibrium energy 

spectrum on eddy viscosity for compressible turbulence 

• Approach: 

^ Use multiple scale concept introduced by 


□ Large-Scale 


_Dk p d w,5L, ,du. 2 

' , sf = a? [( " + ^ ) ^ 1+w( ^ ) p£p + fCl 


p^ = |-[(p+i^)S+ Cp, fpr(§V - Cp 2 p£- + fc 2 
Dt dy (T~ dy k p dy k p 

• fcj - exchanges between the turbulent kinetic energy and internal energy 

• fc 2 - increased spectral energy transfer due to compressibility effects 
□ Small Scale 


_D k t d , HT,dk tl — 

pS = i-[(p+ azffi + - ct,/l 

Dt dy a- dy , k t k t 


□ Eddy Viscosity 


p T r; pul & p[k p + 


i( fc P + fc t) ; 
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Model Evaluation 


• Turbulent Shear Flow 



TURBULENT FREE SHEAR LAYER TURBULENT BOUNDARY LAYER 


• Shock/Turbulent-Boundary-Layer Interactions 
0 transonic flow 



BUMP FLOW BUMP-fN-CHANNEL BUUP-1N-CHANNEL 


0 supersonic flow 



SHOCK REFLECTION RAMP FLOW 


Compressible Turbulent Shear Flow 



TURBULENT FREE SHEAR LAYER TURBULENT BOUNDARY LAYER 
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Bypass transition models 


• Objective: 

0 Develop transition models for flows with free stream 
turbulence 

• Approach: 

^ Using K-e model as the base model 

0 Introduce effective intermittency to either the eddy 
viscosity or the k-e model equations 

Tronsitionol Boundary _cyer 

0.0060 
0,0060 


5 0-0040 

0.0020 


0.0000 

0.0 


0.0000 

0.0060 

D 0.0040 

0.0020 

0.0000 

0.0 





Tronsitionol Boundary Lcyer 
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PDF modeling of turbulent reacting flows 


• Objective: 

^ Develop models that can accurately simulate finite chemical reac- 
tions in turbulent flows. 

0 Develop and validate independent PDF models. 

0 Technology transfer. 

• Approach: 

<£ Joint pdf for scalar compositions. 

<> Moment closure schemes for velocity field. 

^ Develop hybrid solver consisting of Monte Carlo method and 
fini te-difference j finite- volume method. 


BURNER tet 


Without PDF Method 



With PDF Method. (CfflCfTT) 



.M -M .« .11 .13 
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TURBULENCE MODEL DEVELOPMENT AND APPLICATION AT 

LOCKHEED FORT WORTH COMPANY N95- 27884 


Brian R. Smith 
CFD Group 

Lockheed Fort Worth Company 
Fort Worth, Texas 


Broad Range of Flow Problems of Interest 

Wide Range of Flow Conditions: 

Subsonic - Hypersonic 

Internal - External - Store Separation 

Cruise - High Angle of Attack 


Rows phenomena of Interest: 

Inlets/Diffusers 
Streamwlse Curvature 
Shock/BL Interactions 
Rectangular Duct ^ Circular 

Nozzles 

Entrainment 

Round— Rectangular Duct 
High Speed Shear Layers 

External Aerodynamics 
Vortex 

Leading Edge Separation 
Shock/BL Interactions 


Leading Edge Separation - Cowl Lips 
Separation Induced Unstart 


Rim cooling, Liners, Vanes 
Swirl 


3D Boundary Layers 
Wakes 


The CFD Environment at Lockheed Fort Worth Company 


Most codes developed or highly modified in house 
General grid generation and solvers for diverse applications 
Structured and unstructured solvers 
Computational efficiency important 

• Complex geometries, many gridpoints 

• Large arrays of flow conditions 


gBBSJsDtNa 
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Requirements for Turbulence Models 


Turbulence Modeling Priorities for Industrial Application 

• Validation 

High accuracy for attached flows 
Reasonable accuracy for ail flows 
High confidence level 

• Computational efficiency 

• Robust for complex geometries 

• Transitional modeling capability 

To obtain acceptable accuracy, propulsion flows demand more sophisticated 
turbulence models than do external aerodynamic flows 


The k - kl and k - 1 Two Equation Turbulence Models 

Advantages of using kl or I Instead of e or co 

kl and I equations are easier to resolve numerically than c equation 

Dissipation Length Scale is an integral length scale 

•Can derive equation for volume integral of two point correlation function. 
•Theoretical e equation is dominated by small scales 

k - kl and k - 1 agree better with compressible boundary layer data than 
does k~t 

Disadvantage - current formulation requires calculation of distance to walls 


k - kl model k - 1 model 

• Includes unique, consistent wall • Oerived from k - kl model - identical in 

function high Re turbulence 

• Accurate for transonic flows • Near wall model simulates k in viscous 

sublayer 
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The k - kl Model Wall Function 


Wall layer model derived from and consistent with the k - kl model 

• Assume convection In momentum, energy and turbulent kinetic energy 

equations to be negligible 

• Boundary layer approximation 

Match velocity, k and I at first grid point in Navier - Stokes solution 
First grid point can be In viscous sublayer, buffer or logarithmic region 
Boundary conditions on k and I simple for k - kl model 

Advantages of wall functions 

• Reduces number of necessary grid points 

• Reduces number of iterations to converge steady state solution 60 - 90% 


Wall Functions are Accurate for Separated Flow Applications 

Axlsymmetrlc Bump, Transonic Flow Experi ment 


Accurate predictions with and without wall 
functions 
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Velocity profiles with and without wall functions 
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The k - I Model with Near Wall Model 

k! equation is transformed exactly to an I equation 

Advantages of k - I formulation 

• I is linear near wail, kI nonlinear and very small 

• Near wall damping terms disappear 

• Production term drops out with current choice of constants 

k - 1 model includes; 

• Transitional flow modeling 

• Compressibility corrections 

Modeling of details of k profile near wall important for hypersonic flows 

• Magnitude of normal stress term comparable to static pressure 

• Near wall density variations large 


/ Equation Much Easier to Resolve than e Equation 

e equation requires fine grid from wall to y + of 20 to resolve peak 

•Exclusion of near wall viscous dissipation term aggravates problem 
•Logarithmic region, e «= 1 /y 

I equation is nearly linear near wall - much less sensitive to grid resolution 
Length scale and dissipation profiles near waif 



32 


hiiijji 



Resolution Study with k - e and k-l Models 

k - c with Launder - Sharma k -e with Lam - Bremhorst 
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Sample Applications: 


Mach 8 Shock Wave Turbulent Boundary Layer Interactions 


F-16 Inlet Derivitive, Isolated Duct Study 


Multi-slot Ejector 

F1 10 Nozzle Drag Reduction Study 
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k - I Model With Compressibility Correction gives Best Prediction 
For Mach 8 Shock Boundary Layer Interaction 


Mach 8.2 
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Good Predictions of Multi - Slot Ejector Obtained with 

k - kl Model 



Mach Contours 


k - kl Model Predicts Entrainment Effects Near Slots 
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Summary 


Computationally efficient k - 1 and k - kl models have been developed 
and implemented at Lockheed Fort Worth Company 


Many years of experience applying two equation turbulence models to 
complex 3D flows for design and analysis 
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A SUMMARY OF COMPUTATIONAL EXPERIENCE AT GE AIRCRAFT ENGINES FOR 


COMPLEX TURBULENT FLOWS IN GAS TURBINES 

N95- 27885 

R. Zerkle and C. Prakash 
GE Aircraft Engines 
Cincinnati, Ohio 


CONTENTS : 

• INTRODUCTION 

• 2-D BOUNDARY LAYER CODE WITH LRN TURBULENCE MODEL 

• 3-D NAVIER-STOKES CODE WITH WALL FUNCTIONS 

• 3-D NAVIER-STOKES CODE WITH LRN TURBULENCE MODEL 

• FILM COOLING SIMULATION 

• TURBULATED PASSAGE SIMULATION 

• OVERALL CONCLUSIONS 

• LIST OF REFERENCES 


INTRODUCTION : 

• Indications are that the standard k-e turbulence model together with 
standard wall functions are adequate for CFD simulations in cavities 
away from the primary gaspath of a gas turbine engine. 

• However, CFD simulations in the primary gaspath and in blade cooling 
passages require more advanced turbulence models. 

• Therefore, this presentation will summarize some CFD experience at 
GEAE only for flows in the primary gaspath of a gas turbine engine and 
in turbine blade cooling passages. 
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2-D BOUNDARY LAYER CODE WITH LOW REYNOLDS 
NUMBER (LRN) TURBULENCE MODEL : 

• The STAN5 B.L. code was modified to include the LRN k-e turbulence 
model of Lam & Bremhorst as described by Zerkle & Lounsbury [1 ]. 

• Includes the following near-wall effects: 

- High freestream turbulence 

- Axial pressure gradient 

- Onset of transition 

- Relaminarization 

- Wall roughness 

- Wall curvature 

• Used to compute heat transfer coefficient distributions on turbine airfoil 
external surfaces. 

• Primary limitation: 

- It’s a 2-D code in a 3-D environment. 


3-D NAVIE R-STOKES CODE WITH WALL FUNCTIONS : 

• Time- marching finite- volume formulation of the Reynolds -averaged 
Navier-Stokes equations as described by Turner & Jennions [2,3], 

• Includes: 

- Explicit Runga-Kutta flow solver 

- Implicit formulation of the standard k— e turbulence model 

- Standard wall functions 

- Transonic flow effects 

• Used to simulate high speed flows in turbomachinery passages. 

• Limitations: 

- Lacks near-wall physics of the 2-D boundary layer code. 

- For example, lack of boundary layer transition leads to overprediction 
of loss for some turbomachinery airfoil passages containing signifi- 
cant regions of transitional flow. 
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3-D NAV1ER-ST0KES CODE WITH LOW REYNOLDS 
NUMBER (LRN) TURBULENCE MODEL : 

• The LRN k-eturbulence model of Lam & Bremhorst was implementec 
in the 3-D Navier-Stokes code as described by Dailey, Jennions anc 
Orkwis [4]. 

• Addition of the LRN turbulence model improved the prediction of los: 
for transitional flows. 

• Primary limitation: 

- The need for a very fine grid in the near-wall region leads to exces 
sive run times which renders the code impractical for design applies 
tions at this time. 


FILM COOLING SIMULATION : 

• Film cooling at the surface of an HP turbine airfoil is crucial to its life. 

• Improvement of the film cooling process would significantly improve 
turbine performance by reducing the need for cooling air flow. 

• CFD simulation could facilitate film cooling development by reducing 
the need for expensive cascade testing and, more importantly, by giv- 
ing greater insight into the associated flow physics. 

• A CFD simulation of film-cooling tests, which were carried out at the 
Univ. of Texas by Professors Crawford & Bogard, and their students, is 
described by Leylek & Zerkle [5]. 

• These tests are of special interest because the ranges of film cooling 
parameters are consistent with those typically found in gas turbine air- 
foil applications. 

• The objective was to validates CFD model offilmcooling by comparing 
numerical and experimental results. 
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FILM COOLING SIMULATION (CONTD) : 

• The model includes: 

- A 3-D, fully— elliptic, Navier- Stokes solution of the coupled flow in 
the plenum, film hole, and cross-stream regions. 

- An exact representation of the inclined, round, film-hole geometry l 

using a highly- orthogonalized fine grid mesh. J 

- The standard k-e turbulence model with standard wall functions. 

i 





FILM COOLING SIMULATION fCONT’D) : 



Figure 14 . Lateral variation of adiabatic effectiveness from 
computations and experiments for M=0.5 at five 
streamwise stations. 


FILM COOLING SIMULATION (CONTD) : 

• Summary of Results: 

- The flowfield is dominated by a strong three-way coupling between 
the plenum, film-hole, and cross-stream regions. 

- Flow within the film hole is extremely complex, with counter-rotating 
vortices and local jetting effects. 

- A comparison of computed and experimental film effectiveness on 
the plate surface indicates that the simulated coolant jet is not spread- 
ing as fast as experimental results. 

• Conclusions: 

- There is excellent qualitative agreement between the numerical and 
experimental results. 

- However, the lack of lateral spreading of the coolant is caused by the 
inability of the k-e turbulence model to cope with non-uniform rates 
of diffusion in different directions. 

- Improved accuracy requires an anisotropic turbulence model. 
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TURBULATED PASSAGE SIMULATION : 

• Modern high-performance turbine blades are cooled by internal 
radially- rotating serpentine passages. 

• The airflowing through these passages is exposed to very large Coriolis 
and centripetal body forces which induce strong secondary flows and 
buoyant effects. 

• These effects tend to increase heat transfer coefficient on the trailing 
face of an up-pass, but decrease it on the leading face. 

• Turbulators are added to the passage walls in order to enhance their 
cooling effectiveness. 

• The primary objective of blade cooling development is to determine ta- 
bulator and passage configurations which can influence the secondary 
flows to achieve a uniformily high heat transfer coefficient, but within 
pressure-drop constraints. 

• Rotating— passage rig tests are expensive, and it is very difficult to 
achieve high-quality data in the rangeof engine operating parameters. 


TURBULATED PASSAGE SIMULATION (CONT’P) : 

• Therefore, CFD could facilitate blade cooling development by simulat- 
ing new cooling configurations at real engine operating conditions. 

• An exploratory investigation of CFD simulation inturbulated blade cool- 
ing passages is described by Prakash & Zerkle [6]. 

• Conclusions are: 

- The flow fields in turbuiated blade cooling passages are very com- 
plex, and desired accuracy requires advanced turbulence models. 

- An LRN model is needed near turbuiated walls in the case of low pas- 
sage Reynolds number. • 

- An anisotropic turbulent model is needed in the case of large block- 
age ratio (rib height to passage diameter). 

- Practical LRN and anisotropic models are not yet available. 
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OVERALL CONCLUSIONS : 


• Application of the standard k-e turbulence model with wall function 
is not adequate for accurate CFD simulation of aerodynamic perfor 
mance and heat transfer in the primary gas path of a gas turbine engine 

• New models are required in the near- wall region which include mon 
physics than wall functions, the two- layer modeling approach ap 
pears attractive because of its computational economy. 

• In addition, improved CFD simulation of film cooling and turbine blad 
internal cooling passages will require anisotropic turbulence models 

• New turbulence models must be practical in order to have a significar 
impact on the engine design process. 

• A coordinated turbulence modeling effort between NASA center 
would be beneficial to the gas turbine industry. 
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THE APPLICABILITY OF TURBULENCE MODELS TO AERODYNAMIC AND PROPULSION 
FLOWFTELDS AT McDONNELL DOUGLAS AEROSPACE 

Linda D. Krai, John A. Ladd, and Mori Mam 
McDonnell Douglas Aerospace 
St. Louis, Missouri 



Objective 

• Evaluate turbulence models for integrated aircraft components 
such as the forebody, wing, inlet, diffuser, nozzle, and afterbody 

Approach 

• Integrate turbulence models into existing Navier— Stokes program 
maintaining zonal philosophy 

• Introduce corrections to baseline turbulence models to account 
for additional affects such as compressibility or separation 

• Develop algorithmic improvements for better numerical stability 
and robustness 

• Compare the strengths and weaknesses of turbulence models 

• Determine applicability of algebraic, one-equation, and two-equation 
turbulence models for typical complex flows and geometries 
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Turbulence Modeling Capabilities 

• Algebraic Models 

- Cebeci-Smith boundary layer model 

- Baldwin-Lomax boundary layer model 

- P. D. Thomas shear layer model 

• One-Equation Models 

- Baldwin-Barth 

- Spalart-Allmaras 

• Two-Equation Models 

- High Reynolds number k — e 

- Low Reynolds number k — e (Jones-Launder, Speziale, Chien, 
Lam-Bremhorst, So, and Huang-Coakley) 

- Wilcox k — lu 

- Menter baseline and shear-stress transport blended k — w/k — e 


Navier-Stokes Time-Dependent Algorithm 

NASTD 


• Euler /Navier-Stokes Equations 

- Laminar or Turbulent 

- Ideal Gas, Thermally Perfect Air, Equilibrium or Nonequiiibrium Chemistry 

• Finite Volume Formulation 

- Roe and Coakley Flux Difference Split Schemes, Optional TVD Schemes 

• Solution Update Procedure 

- Approximate Factorization 

- Runge-Kutta Time Stepping 

- Iterative Space Marching (PNS) 

• Geometric Capabilities/Generalizations 

- Zonal Capabilities and Flexible Boundary Conditions 

- Grid Sequencing 

- Overlapping Grids 

• Turbulence Models 

- Cebeci-Smith, Baldwin-Lomax and P. D. Thomas Algebraic Models 

- Baldwin-Barth and Spalart Allmaras One-Equation Models 

- Six Low Reynolds Number k - e Models 

- k - u; and Menter blended k - uj/k — t Models 
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Selected Applications 


• Transonic Supercritical Airfoil 

• Three-Element High-Lift System 

• Single Slot 2-D Ejector Nozzle 

• Confluent Mixer 

• Highly Offset 3-D Diffuser 


Modifications to Production Term 


Default calculation of production: 

1 Idui <9tiA 2 2/dufc\ 2 

2\dxj dx{ ) 3\<9 xj 

Vorticity used in production: 




2 _ du k 

r k a^ 


Production limiter used: 

= min(Pjt, 20Dfc) = min(Pfc, 20 c 2 p k Re) 
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Effect of Production Limiter for the Chien k-e Model 

RAE Airfoil Analysis, Turbulent Viscosity Contours 
Mach = 0.725, a = 2.55 deg., Re = 6.5 Million 



RAE Airfoil Analysis 

Moo = 0.725, a = 2.55°, Re = 6.5 Million 

Effect of Freestream Turbulence Level on Surface Pressure 
Cbien k — e Turbulence Model 



X 
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RAE Airfoil Analysis 

Moo = 0.725, a = 2.55°, Re = 6.5 Million 

Production Limiter Used 
Chien k — e Turbulence Model 



RAE Airfoil Analysis, Turbulent Viscosity Contours 
Mach = 0.725, a =2.55 deg., Re = 6.5 Million 
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RAE Airfoil Analysis, Mach Contours 
Mach = 0.725, a =2.55 deg., Re = 6.5 Million 



IIIIHI Mil ,11 1 till i |i| I || |,|| , HI (I), lillinl 



NASTD Solutionof MDA Three-Element High-Lift System 


M-OA AOA-16J21 














Four- Zone Grid for an Ejector Nozzle 



ir-m-L. i. t 3=i -T t==hd 

T T 

Primary Nozzle Symmetry Plane 


Single Slot Ejector Analysis 
NPR=14., Pts/Ptp=.34 

Mach Number Contours from Several Turbulence Models 
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Single Slot Ejector Analysis 

NPK=14., Pts/Ptp=.34 

Eddy Viscosity from Several Turbulence Models 




Single-Slot Ejector Nozzle Analysis 

NPR = 14 , P t JP t = 0.34 

h!h - 

Comparison of Predicted Ejector Flow Rates 


Model 

w,/w p 

% Error 

Experiment 

0.1010 

— 

Thomas /Baldwin-Lomax 

0.1108 

+9.7 

Baldwin-Barth 

0.1129 

+ 11.8 

Spalart-Allmaras 

0.1146 

+13.5 

i 

Chien k — e 

0.1168 

+ 15.6 

Jones-Launder k — € 

0.1126 

+ 11.5 

Speziale k — e 

0.1127 

+11.6 

So k — e 

0.1148 

+13.7 

Huang-Coakley k — e 

0.1112 

+ 10.1 
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Single Slot Ejector Nozzle 

Surface Static Pressure Comparison with Experimental Data 
NPR = 14.0, PJP lp = 0.34 




Offset Diffuser Analysis 

Ae/At=1.6, L/D=4.5, Design Pressure Ratio 

Surface Pressure and Computational Mesh 
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Offset Diffuser Analysis 

Ae/At=1.6, L/D=4.5, Design Pressure Ratio 
Comparison of Engine Face Total Pressures 




Baldwin-Lomax 


Baldwin-Barth 


Spalart-Allmaras 




Offset Diffuser Analysis 

Ae/At=1.6, L/D=4.5, Design Pressure Ratio 
Comparison of Engine Face Total Pressures 



Experiment ! Ebien 
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Offset Diffuser Analysis 

Lower Centerline Surface Static Pressure 
A e /A t = 1.6 , L/D = 4.5, Design Pressure Ratio 



x (cm) 


Offset Diffuser Analysis 

Upper Centerline Surface Static Pressure 
A e /A t = 1.6, L/D = 4.5, Design Pressure Ratio 
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Three-Dimensional Highly Offset Diffuser 

A e jA t — 1.6, LjD — 4.5, Design Pressure Ratio 
Comparison of Engine Face Parameters 


Model 

Pt CV3 / Pi <«, 

P t mi JPtoo 

P Uvo 

Experiment 

0.958 

0.890 

0.114 

Baldwin-Lomax 

0.936 

0.708 

0.292 

Baldwin-Barth 

0.944 

0.735 

0.265 

Spalart-Allmaras 

0.955 

0.860 

0.140 

Chien k — e 

0.970 

0.894 

0.106 

Jones-Launder k - e 

0.966 

0.896 ! 

0.104 ! 

i 

So k — e 

0.975 

0.888 

l 

0.112 
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Comparison of Throat Total Temperatures 

BL=2.5 BL=1.5 BL=0.5 BL=0.5 BL=1.5 BL=2 



Experiment Spalart-Allamaras 


GE Slot Cooled Nozzle, Confluent Mixer 

Surface Temperature Distributions, 

TGX = (T, - - T, m ) 



Baldwin- Barth Model 

k-Epsiion Model 

Spalart-Allamaras 

O Experiment (Average) 
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GE Slot Cooled Nozzle, Confluent Mixer 

Surface Temperature Distributions, 

TGX = (Tr - T Jcoli )/(T Ihol - T Ico j 





Summary of Turbulence Modeling 
at McDonnell Dougalas Aerospace 

• The one-equation models have replaced the algebraic models as the 
baseline turbulence models. 

• The Spalart-Allmaras one-equation model consistently performs bet- 
ter than the Baldwin-Barth model, particularly in the log-layer and 
free shear layers. Also, the Spalart-Allmaras model in not grid de- 
pendent like the Baldwin-Barth model. 

• No general turbulence model exists for all engineering applications. 

• The Spalart-Allmaras one-equation model and the Chien k-e mod- 
els are the preferred turbulence models. 

• Although the two-equation models often better predict the flowfield, 
they may take from two to five times the CPU time. 

• Future directions are in further benchmarking the Menter blended 
k - ujk - e and algorithmic improvements to reduce CPU time of 
two-equation model. 
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EXPE RIEN CE WITH k-e TURBULENCE MODELS FOR HEAT TRANSFER 

COMPUTATIONS IN ROTATING N95- 27887 


Prabhat Tekriwal 

GE Corporate Research and Development 
Schenectady, New York 


OUTLINE 

• Geometry and flow configuration 

• Effect of y+ on heat transfer computations 

• Standard and Extended k-e turbulence model results with wall 
function 

• Low-Re model results (the Lam-Bremhorst model without wall 
function) 

• A criterion for flow reversal in a radially rotating square duct 

• Summary 



Fig. 1-lllustration of geometry anti physics of flow 
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TWO-EQUATION TURBULENCE MODELS 


Mt = f(t q, p k% 


D(pk) _ J)_ ll|_ (Ik 
Dl ” 3xj 


Dx'i» + P (r * ' e) 


DfeD 

Dl 


3 Ht 3e e Gt 

Si <f k £? +f|C 'k p c * - & c 2 p ^ ♦ c 3 p 


where 


_ [it ,duj 3ui 3uj „ 

G k= p <sj + si>sj- c i‘= 0 - 09 


Standard k-c model: 

Prk = 1.0, Pr e =1.3, Ci=1.44, 02=1.92, C3=0.0,fi=1.0. f 2 = 1 . 0 , and fn=1.0 
Extended k-e model: 

Prk = 0.89, Pr e =1.15, Ci =1.1 5, C 2 =1.9, C 3 =0.25 , f,=1.0, f 2 =i.O, and f,,=1.0 


I 


Lam-Bremhorst low-Re model: 

Prk =1.0, Pr e =1.3, Ci =1.44, C 2 =1.92, C3=0.0, f t =(1 +0.05/f|,) 3 , f 2 =1-e R t, 
and fii=(1-e'0 0165R k ) 2 (i +20.5/RO, where Rk=k 1/2 y p/ji and R t =k 2 p/p e 
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Nu/Nu°° Nu/Nuoo 



Fig. 4(b)-Comparlson of model results with data on trailing wall 



Fig. 4(c)-Comparlson of model results with data on leading wall 
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Fig. 6(a)-Plots of Nu-ratlo and k-ratlo from the two models 
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Fig. 7 (a)-Comparison on both leading and trailing walls 



Fig. 7 (b)-Comparison on both leading and trailing walls 


73 




Fig. 8(a)-Comparlson of Eke and low-Re results on leading wall 
(Re=1 0000, Ro=.088, Case B) 



Fig. 8(b)-Comparlson of Eke and low-Re results on leading wall 
(Re=500D,Ro=.1 76,Ca$eB) 
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z/d 

Fig. 8(c)-Comparlson of Eke and low-Re results on leading wall 
(Re=5000,Ro=.176,CaseC) 



Fig. 8(d)-Comparlson of low-Re results with data (Re=2500,Ro=.352,CaseB) 
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Table I Prediction of Row Reversal Near Ihe Leading Wall 


Ro 

AT/Tw 

R/d 

Rc 

Gr/Rc* 

Flow Reversal ? 

0.12 

0.07 

49 

25000 

0.05 

mm 





0.09 

■ ■ 





0.16 

■ 


0.36 



0.26 



0.48 



0.34 

19 


0.07 

196 


0.20 

No 



300 


0.30 

Yes 

0.24 

0.07 

33 

25000 

0.13 

No 



49 


0.20 

No 



196 


0.77 

Yes 



300 


1.18 

Yes 


0.13 

49 


0.36 

Yes 


0.16 



0.45 

Yes 


0.23 



0.65 

Yes 


0.07 


12500 

0.20 

No 


0.13 



0.36 

Yes 


0.23 



0.65 

Yes 

0.34 

0.13 

49 

25000 

HEIQHi 

Yes 


0.16 




Yes 


0.23 



8» MB 

Yes 

0.48 

0.13 

49 

25000 

1.45 

Yes 


SUMMARY 


1. Near-wall grid size has a significant effect on the heat transfer calculations 
when the "wall function" treatment is used. Numerical experiment on the data 
of Morris et al. (1991) suggests that a y+ value in the range of 1 2 to 42 or so 
yields more accurate results. 

2. The extended k-e turbulence model, while yielding heat transfer results virtually 
the same as those of standard k-e model for low rotation-number flows, 
provides an improvement over the standard k-e model by up to 1 5% or so in 
heat transfer predictions for high rotation number flows. 

3. Wall-function k-e models predict lower (than data) heat transfer at the trailing 
wall and higher at the leading wail. The need to properly represent the effect of 
rotation in the k-e model equations is realized. 

4. The low-Reynolds number model utilizes a large number of cells and the 
convergence rate is very slow in comparison to the high-Reynolds number 
model using wall function. It is difficult and expensive to obtain a well 
converged solution with the low-Re turbulence model. 
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5. The poor agreement of the low-Re model results with the data makes the low- 
Re model as an unattractive choice for heat transfer computations in rotating 
radial outward flow at high Rotation number (> 0.24) and high-Reynotds 
number (25000). 

6. The extended version of high-Reynolds number turbulence model in 
conjunction with wall function yields satisfactory results for flows with 
isothermal walls as well as uneven wall temperatures. The agreement is within 
5-25% of the data with uneven wall temperatures for flows at Reynolds 
numbers 10000 or higher. 

7. For flows at Reynolds number 5000 or lower, the low-Re model predictions are 
better, especially for the case of uneven wall temperature conditions. 

8. The centrifugal buoyancy may cause a flow reversal near the leading wall 
depending upon the geometry and flow parameters such as rotation number, 
temperature ratio, mean radius ratio and Reynolds number. For the square- 
section channel considered here, a criterion of Bo=Gr/Re2 higher than 0.3 is 
predicted to cause flow reversal near the leading wall for flows at Reynolds 
number up to 25000. 
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CONTENTS 


• Gas Turbine Combustor Flow Physics 

• Turbulence Model Investigations 

• Turbulent Combustion Modeling 

• Present Status and Future Needs 


GT COMBUSTOR FLOW PHYSICS 


• Key issue is flame stabilization by means of recirculating flow 
of hot gases and chemically-active species to ensure continuous 
ignition of fresh reactants. 

• Three main mechanisms: 1) axial swirling air jet associated with 
each fuel introduction; 2) sudden expansion of axial swirling 
jets; 3) blockage due to radial air jets downstream of fuel sources. 
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TURBULENCE MODELS SURVEYED 


• Following models or modifications have been tested at P&W / 
UTRC using RANS solvers on building block flows* 

1. low- Re models (complex ducts); 

2. RSTM or SMC (complex ducts, swirling and non-swirling 
dump combustor); 

3. RNG (pipe, backstep, 180 deg duct); 

4. two-layer near-wall model (internal flows, heat transfer); 

5. realizable algebraic stress model (swirling dump combustor); 

6. compressible turbulence (shear layers, compression corner) 

7. steady vs. unsteady-state solver (bluff-body, compression 
corner) 

• Major difficulty occurs with swirling flows, and failure to predict 
downstream velocity components. 


SWIRLING FLOWS 


• Benchmark-quality data set provided by Johnson-Roback 
co-annular combustor with swirl: 




• Poor agreement of CFD and data highlights need for improved 
upstream BC specification (swirler geometry), 3-D, unsteady 
analysis. Even SMC models fail to reproduce downstream 
velocity profiles. 
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UNSTEADINESS AND FLOW FIELD RESOLUTION 

• RANS solvers can predict flow coherence (vortex shedding) 
when run in an unsteady mode with small At. 

• Same flow field computed in steady-state sense gives completely 
unusable results. 

• Example: V-gutter flow, computed by Durbin (1994): 


U-COMTOURS 

STEADY-STATE COMPUTATION 



3.4 


UNSTEADY -STATE COMPUT ATION 



UNSTEADINESS AND FLOW FIELD RESOLUTION 

• RANS solvers cannot predict flow oscillations at frequencies 
near characteristic turbulence frequency. 

• Example: Unsteady comp, comer flow of Dolling and Or (1983): 

Frequency ratio: 


M. - 2M T t - 486' R 
& - LQ2 in. Re 5 - 1.78 X 10* 

• Separation bubble oscillations (at resonant frequency) not 
resolved by RANS solver. 

• Limitations of steady-state and unsteady-state RANS solvers set 
by flow characteristic time scales. 

True time-accurate solvers (LES, DNS) needed for prediction 
of all relevant phenomena 
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TURBULENT COMBUSTION MODELING 

• Eddy Dissipation Concept Model, together with reaction exclusion regions, capable 
of prediction gross flow features at near LBO conditions (Sturgess et al., 94-GT-433) 



Reaction exclusion regions 



AjdMOhUnea 


EDC-mode] Temperature field 


• EDC model, however, fails to predict flame attachment at rich conditions 


TURBULENT COMBUSTION MODELING 

• Assumed-Pdf method of Girimaji (LaRC Workshop, 1991) used with 
non-equilibrium kinetics model. 



Temp, °R — * Temp, °R -*■ 


• Example: N + O 2 ^5 NO + 0 in extended Zeldovich model 

• Results dependent on T Low , T ffigh , 4>, modeling of hh transport equation, etc. 

• More testing needed 
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PRESENT STATUS OF COMBUSTOR MODELING 

• Corsair (Ryder, P&W) unstructured, unsteady flow solver 



• Example: Time-dependent combustor flow using engineering 
boundary conditions, compressor exit to turbine inlet 

• Code currently includes standard k-e and EBU combustion 
model. Additional capabilities being added under ’’Subsonic 
Emissions and Combustor Design Code” program with NASA LeRC. 


PRESENT STATUS OF COMBUSTOR MODELING 


• Example: Structured flow solver solution of Task 200 
LBO Research Combustor: 



Temperature 


300.0 


1275.0 


2250.0 


3225.0 


4200.0 


• k-e turbulence model 

• EBU combustion model for propane fuel 

• 285,000 elements 
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PRESENT STATUS OF COMBUSTOR MODELING 


• Example: Unstructured flow solver solution of Task 200 
LBO Research Combustor: 



• EBU combustion model for propane fuel 

• Approx. 300,000 elements 


TURBULENCE RESEARCH NEEDS 


• Modelling: Applications / validations of currently available 
combustion models (P*pdf, Monte Carlo pdf, laminar flamelet) 
to complex combustor geometry with jet fuel kinetics. 

• Flow Physics: Accurate numerical description of mechanisms 
responsible for flame holding, local extinction (LES, DNS); 
contrast cold flows with heat release flows. 



Entrainment of unbumed fuel 
in the recirculation region 
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COMBUSTION SYSTEM CFD MODELING AT GE AIRCRAFT ENGINE S 

N95- 27889 

D. Burrus and H. Mongia 
GE Aircraft Engines 
Cincinnati, Ohio 


and 


A. Tolpadi, S. Correa, and M. Braaten 
GE Corporate Research and Development 
Schenectady, New York 


CURRENT COMBUSTION SYSTEM CFD MODELING CAPABILITIES AT GEAE 
PROVIDED BY THE CONCERT CODE 

KEY FEATURES INCLUDE; 

FINITE VOLUME, PRESSURE CORRECTION FORMULATION 
SECOND ORDER ACCURATE QUICK NUMERICS 
SINGLE STRUCTURED BODYFITTED GRID 

CONVENTIONAL K-E TURBULENCE MODEL WITH LOG WALL FUNCTIONS 

AVAILABLE COMBUSTION MODELS INCLUDE; 

SINGLE SCALAR PRESUMED SHAPE PDF (FAST CHEMISTRY) 

TWO SCALAR PRESUMED SHAPE PDF (REACTION PROGRESS VARIABLE) 

TWO STEP EDDY BREAKUP (ARRHENIUS KINETICS) 

ZELDOVICH THERMAL NOx MECHANISM (FORWARD AND REVERSE REACTIONS) 

BOTH 2D/AXISYMMETRIC AND FULLY 3D VERSIONS AVAILABLE AND IN DAY TO DAY 
USE 

CURRENTLY HAVE A USER BASE OF OVER 20 ENGINEERS AT GEAE AND GE-CRD 

TYPICALLY APPLIED TO PREDICT COMBUSTOR PERFORMANCE INCLUDING; 
EMISSIONS (CO, HC, AND THERMAL NOx), COMBUSTION EFFICIENCY 
EXIT GAS TEMPERATURE RADIAL PROFILE AND PATTERN 
GENERAL FLOW FIELD CHARACTERISTICS 
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CONCERT DEVELOPMENT HISTORY 


EFFORT INITIATED IN 1983 

INITIAL PRODUCTION VERSION RELEASED TO GEAE USERS IN 1987 

FOCUSED TO PROVIDE HIGHLY PRODUCTIVE ENGINEERING ANALYSIS CAPABILITIES 

- GRID GENERATION OPTIMIZED FOR THE SPECIFIC GEOMETRY FEATURES OF THE GAS TURBINE 
COMBUSTOR 

• - INCLUDES ROUND DILUTION HOLES, SWIRLER DISCHARGE, AND LINER SLOT FEATURES 
WITHIN THE GRID 

- EASY INTRODUCTION OF INTERNAL BODIES OF COMPLEX GEOMETRY 

- WORKSTATION BASED USER FRIENDLY PRE AND POST PROCESSING FUNCTIONS BUILT AROUND 
THE SOLVER 

- SOLVER HIGHLY OPTIMIZED FOR THE GEAE CRAY C-90 COMPUTER 


TYPICAL ID MODEL OF A COMBUSTOR UTILIZING A MESH OF -100,000 POINTS CAN BE GENERATED, RUN, AND POST 
PROCESSED WITHIN A SINGLE WORKING DAY I 


HAS UNDERGONE CONTINUAL DEVELOPMENT TO IMPROVE AND ENHANCE MODELING 
CAPABILITIES 


| 


- CURRENTLY ON VERSION 3 RELEASE 
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COMBUSTION SYSTEM CFD MODELING IN ACTION AT GEAE 



SWIRLCUP/SPRAY MODELING 

• RECIRCULATION SIRENGTH/SIZE 

• FLOW FIELD CHARACTERISTICS 

• STRAY DROPLET TRAJECTORIES 

• INPUTS FOR 3D COMBUSTOR MODEL 



DIFFUSER FLOW MODELING 

• P* RECOVERIES AND Pt LOSSES 

* FLOW FIELD CHARACTERISTICS 





MODELING APPLIED FOR DESIGNING ENGINE COMBUSTION SYSTEMS 


PRODUCTION ENGINES DEMONSTRATOR ENGINES 


CFM56-5B DUAL ANNULAR 
GE90 

CF6-80C LOW EMISSIONS 
LM1600 DLE 
LM2500DLE 
LM6000 DLE 


YF120 

FI20 

XTE45 IHPTET PHASE I DEMO 
XTE46 IHPTET PHASE H DEMO 


ADVANCED ENGINES 


A/F-X 

NASA/GE HSCT 

NASA ASI PRELIMARY CONCEPTS 
DOE/GE ATS 


MODELING APPLIED TO IMPROVE FUNDAMENTAL UNDERSTANDING 


CFM56-3 AND CFM56-5B NOx EMISSIONS CHARACTERISTICS DIFFERENCES 

CFM56-5A EXIT GAS TEMPERATURE PROFILE SHIFT 

FI20 PATTERN FACTOR AND RADIAL PROFILE IMPROVEMENT 

LM2500 CO EMISSIONS REDUCTION EFFORT 

CF34 LINER COOLING MOD IMPACT ON CO EMISSIONS 

FUOX AUGMENTOR MIXER, SPRAYBAR, FLAMEHOLDER INTERACTION OPTIMIZATION 
F1KW00 AUGMENTOR EXHAUST DUCT LINER FAILURE AND FIX INVESTIGATION 


gNOx/kgfciel 







NON-DIMENSIONAL EXIT PASSAGE HEIGHT 


CONCERT3D vs. RIG DATA COMPARISON FOR NASA/GE E3 COMBUSTOR 
' (EXIT GAS TEMPERATURE AVERAGED AND MAXIMUM RADIAL PROFILES) 



CONCERT3D vs. RIG DATA COMPARISON FOR NASA/GE E3 COMBUSTOR 

(NOx EMISSIONS) 
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GEAE CONCERT EXPERIENCE: 


CONCERTED WITH PRESUMED SHAPE PDF/FAST CHEMISTRY MODEL AND 
THERMAL NOx MODEL DOES WELL AGAINST REAL ENGINE DATA 


CONCERT3D WITH TWO STEP EDDY BREAKUP MODEL DOK NOT ^ 
CONSISTENTLY DEMONSTRATE ACCEPTABLE AGREEMENT FOR [CO] 

AND ^ HC ^ ( ^^ S pg^ F0RMANCE ISSUES NOT AS WELL PREDICTED COMPARED 
TO PRESUMED SHAPE PDF/FAST CHEMISTRY APPROACH 


SHORTCOMINGS: 

TWO STEP EDDY BREAKUP MODEL NOT ADEQUATE FOR THE REQUIRED 
LEVEL OF PREDICTIVE ACCURACY 


FAST CHEMISTRY CANNOT PREDICT [CO], [HC], AND IGNITION, 
BLOWOUT, AND RELIGHT 


REQUIRES ACCURATE FINITE RATE CHEMISTRY REPRESENTATION 
AND MORE ACCURATE TURBULENCE-CHEMISTRY INTERACTION 


MODELING 


CE HAS EMBARKED ON THE DEVELOPMENT OF IMPROVED CONCERT MODELING CAPABILITIES 


HYBRID CONCERT CFD / MONTE-CARLO MODELING APPROACH 

APPROACH ADOPTED FOR THE NEXT RELEASE OF COMBUSTION CFD MODELING 
CAPABILITY AT GEAE 


RETAINS; 


- SINGLE STRUCTURED BODYFITTED GRID 

- PRESSURE CORRECTION FINITE VOLUME FORMULATION 

- K-E TURBULENCE MODELING WITH LOG WALL FUNCTIONS 


INTRODUCES; 

- MONTE-CARLO SCALAR PDF TO ADDRESS TURBULENT COMBUSTION 


- SINGLE ATTRIBUTE (CONSERVED SCALAR) FOR FAST CHEMISTRY 

- MULTIPLE SCALARS FOR FINITE RATE CHEMISTRY OF CH4 AND JETA FUELS 
BASED ON APPROPIATE REDUCED MECHANISMS 


DEVELOPMENT HAS BEEN UNDERWAY SINCE 1992 


- 3D CODE DEVELOPMENT INITIATED IN' MID YEAR 1993 
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HYBRID CONCERT CFD / MONTE-CARLO MODELING APPROACH 



SCHEMATIC OF COMMUNICATIONS IN THE COMBINED CONCERT / MONTE-CARLO MODELING 


HYBRID CONCERT CFD / MONTE-CARLO MODELING APPROACH 

BETA TESTING INITIATED BEGINNING OF 1994 

FOCUSED ON FAST CHEMISTRY CALCULATIONS AND OPTIMIZING 
COMPUTATIONAL EFFICIENCY 


SIGNIFICANT IMPROVEMENT IN COMPUTATIONAL EFFICIENCY ACHIEVED 

TEST CASE 1 TEST CASE 2 


NUMBER OF GRID POINTS 
NUMBER OF M/C PARTICLES 
CPU TIME (CRAY C-90 seconds) 

9,261 

216,000 

58,621 

1,500,000 

CONCERT WITHOUT M/C 

83 

5,400 

INITIAL HYBRID CONCERT /MC 
OPTIMIZED VERSION 
PERCENT REDUCTION 

39,960 

1,770 

-95.6% 

187,560 

41,400 

-77.9% 

WALL CLOCK TIMES (seconds) UTILIZING CRAY 
MULTI-TASKING OPTION 

1,500 

29,520 


RUN TIMES HAVE BEEN REDUCED TO THE POINT WHERE OVERNIGHT TURNAROUND TIMES FOR A TYPICAL 
3D COMBUSTOR MODEL ARE POSSIBLE 
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HYBRID CONCERT CFD / MONTE-CARLO MODELING APPROACH 

(INITIAL 3D CALCULATION OF CFM56-3 COMBUSTOR WITH FAST CHEMISTRY) 



CALCULATED FLOW FIELD IN PLANE IN LINE WITH INLET SWIRL CUPS 

INITIAL CALCULATED RESULTS SHOW A TEMPERATURE FIELD THAT DOES NOT AGREE WELL WITH EXPECTED 
LEVELS. CALCULATION SHOWS CONSIDERABLY LESS DIFFUSION OF THE SCALAR FIELD (FUEL MIXTURE 
FRACTION) THAN OBSERVED FROM RIG DATA AND CONCERT CALCULATIONS PERFORMED USING THE PRESUMED 
SHAPE SCALAR PDF COMBUSTION MODELING APPROACH. 


HYBRID CONCERT CFD 7 MONTE-CARLO MODELING APPROACH 


FUTURE WORK PLANNED 

- PERFORM CALCULATIONS AGAINST A BENCHMARK REACTING FLOW 
EXPERIMENT WITH AVAILABLE TEST DATA 

- BLUFF BODY STABILIZED FLAME ; (GULATI AND CORREA) 

- SYSTEMATICALLY STUDY THE EFFECTS OF SCHMIDT NUMBER AND OTHER 
PARTICLE TRACKING PARAMETERS ON THE FAST CHEMISTRY SOLUTION TO 
IMPROVE AGREEMENT WITH THE DATA 

- PERFORM 3D SINGLE AND DUAL ANNULAR COMBUSTOR CALCULATIONS AND 
COMPARE RESULTS WITH AVAILABLE GEAE DATA BASE 

- IMPLEMENT REDUCED CHEMISTRY SCHEMES (MULTIPLE SCALARS) TO PERFORM 
FINITE RATE CHEMISTRY CALCULATIONS 

- PREDICT (CO], [HC], AND [NOx] EMISSIONS 

- RELEASE CODE FOR PRODUCTION USE AT GEAE 

- FAST CHEMISTRY BY END OF FIRST QUARTER OF 1995 

- FINITE RATE CHEMISTRY BY END OF THIRD QUARTER OF 1995 
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FUTURE MODELING DIRECTIONS 


FOCUSED ON IMPROVING THE PREDICTIVE ACCURACY FOR ALL KEY 
COMBUSTOR PERFORMANCE ISSUES TO LEVELS THAT WOULD ELIMINATE 
THE NEED FOR COMPONENT RIG DEVELOPMENT TESTING 


1970’s / 1980’s 



WEEKS ► 


FUTURE MODELING DIRECTIONS 


INDUSTRY WILL LOOK INCREASINGLY TO THE ACADEMIC COMMUNITY 
(UNIVERSITIES AND NATIONAL LABS) TO DEVELOP THE NEEDED 
MODELING IMPROVEMENTS 


INDUSTRY MUST PROVIDE THE GUIDANCE AS TO WHAT IS NEEDED 


FUTURE GENERATION MODELS MUST) 


- PROVIDE MORE RIGOROUS REPRESENTATION OF COMPLEX PHYSICAL 

- BE CCOTEFEECTIVE AS A ROUTINE APPLIED DESIGN/ANALYSIS TOOL 

I ^vroETOE^^OFAC^S^^APABILrr^ DEMANDED 
OF IT 


COMPUTING PLATFORM CAPABILITIES ARE ADVANCING AT A RAPID 
PACE 


THE PRACTICALITY OF ADVANCED MODELS IN INDUSTRY MAY NOT BE 
TOO FAR INTO THE FUTURE 


TIMF TO START NOWON DEVELOPMENT OF THE ADVANCED MODELS OF THE FUTURE INTO PRACTICALTOOLS 
TO HAVE THEM READY FOR USE WHEN THE REQUIRED COMPUTING PLATFORMS BECOME AVAILABLE IN 
INDUSTRY 
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Calculation of turbulent heat transfer 
in "cluttered spaces", by BRIAN SPALDING 
Topic is The WDIS & WGAP calculation. 


The need: 

* Prandtl-raixing-length models require knowledge of distance 
from nearby walls AND between walls (eg Nikuradze formula) 

* Many low-Re models require the distance from nearby walls 

* In spaces "cluttered" with solids (eg electronics cooling) , 
calculation of distances and gaps has, in the past, been 
time-consuming. 

The solution: 

* This contribution computes WDIS and WGAP (the required 
quantitities) by solving: 

div grad L = -1 

with L fixed to zero in solids. 


— J 

Computational 

Turbulence 

1994 


Outline of the theory 


Obviously L values which satisfy this equation will be proport- 
ional to the distance from the wall at points which are close 
to it. The question is: what is the proportionality constant? 

The constant depends also on the distance across the inter- 
solid space, which however is the other unknown which it is 
desired to determine. 


The practice adopted by the author is to deduce both the 
required quantities, WDIS the distance from the wall, and 
WGAP the distance between walls (whatever these quantities may 
mean in "cluttered spaces"), from the an algebraic fucntion 
of the local values of L and its gradient. 


— — 1 
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1994 

3 
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The results 


The formula employed gives exact results for situations where 
WDIS and WGAP have unequivocal meanings, namely for the space 
between two parallel plates or within a long circular-sectioned 
pipe; and it gives plausible results for more complex cases. 

The equation for L, with the appropriate boundary conditions, 
is of course very easy to solve by numerical means; so WDIS 
and WGAP can be quickly computed before the flow simulation 
starts. 

The use of the method is illustrated by a PHOENICS calculation 
for a geometry involving two boxes, a connecting arc, an inlet 
and an outlet. It was performed by I Poliakov and S Semin, 
of CHAM, to whom the author's thanks are due. 
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1994 


Calculation of turbulent heat transfer 
in "cluttered spaces", by BRIAN SPALDING 
Topic 2 . The LVEL model. 


The need: 

* In "cluttered" regions, the between-solid distances are 
too often too small for fine-grid resolution. 

* Reynolds numbers are usually low, at least in some plsces. 

* A model is needed which gives plausible results in these 
circumstances AND fits experimental data for better-studied 
ones. 

The solution: 

* The LVEL model of PHOENICS gets local effective viscosities 
from the analytical nuplus -versus -up lus relation which fits 
the laminar, transitional & full-turbulent ranges very well 
Only local velocity and WDIS (wall distance) are needed. 


1 1 
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Outline of the theory 


The u-plus versus y-plus formula of Spalding (1961) is employed 
namely: 

y+ — u+ + { 1/E) * [ exp{K*u-f ) - 1 - K*u+ - (K*u+)**2/2 

- (K*U+)**3/6 - (K*u+)**4/24 ] 

which implies the formula foi: dimensionless effective viscosity 
V+ =14- (K/E) * ( exp (K*U+) - 1 - K*U+ - (K*U+)**2/2 

- (K*U+)**3/6 ] 

With the wall-distance and the velocity known at every point, 
the effective viscosity can also be computed at every point. 

The method is valid for the whole range of Reynolds numbers; 
but it is best supplemented by a low-Re ,f v+-collapse M formula. 


Computational 

Turbulence 

1994 



The results 


The LVEL model gives the well-known experimental results for 
simple circumstances, such as flow in pipes and between parall- 
plates; and it gives plausible results for more complex cases. 

The use of the method is illustrated by a PHOENICS calculation 
of the flow and heat transfer in a small part of a large space 
cluttered with solids which participate in the heat-transfer 
process . 

The method is the only plausible and practicable one known to 
the author for handling heat transfer in electronics-cooling 
problems, because of the excessive grid-fineness requirements 
of low-Reynolds-number k-epsilon extensions. 
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RECENT PROGRESS IN THE JOINT VELOCITY-SCALAR PDF METHOD 


M.S. Anand 

Allison Engine Company 
Indianapolis, Indiana 


o TURBULENCE 

o REACTION (treatment, kinetic schemes, emissions) 
o TURBULENCE/CHEMISTRY INTERACTIONS 
o ATOMIZATION 
o SPRAY EVAPORATION 

SIMULATION ISSUES: 
o NUMERICS (accuracy, convergence) 
o GEOMETRY (body-fitted grids, unstructured grids) 
o COMPUTATIONAL RESOURCES (Time, Storage) 
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JOINT VELOCITY-SCALAR PDF METHOD 


SIGNIFICANT MILESTONES AND RECENT PROGRESS 

o 2-D and 3-D time dependent flows (with finite-volume method) 

(Anand et al. 1987, Haworth & El Tahry 1989) 

o Stochastic dissipation model development and validation 
(Pope & Chen 1990, Pope 1991, Anand et al. 1993) 

o 2-D Elliptic flows (mean pressure algorithm), swirling flows 
(Anand et. 1989, 1993) 

o Spray treatment 
(Anand 1 990) 

o Manifold methods for reaction kinetics 

(Maas & Pope 1992, 1994; Norris & Pope 1994; Norris & Hsu 1994) 


o Solve Poisson equation for mean pressure: 


9 z <p> 
j j 


a 2 

dx.dx. 

> J 


<pU i U j > 


o Satisfy continuity by solving for velocity correction potential, velocity correction 

1 3<j> 

dx. 


3 2 <|> 


dxdx. 


d TT 

■'5T <pU - > 


AU. = ■ 


<p> 


o Solution algorithm is consistent with B-spline representation of mean fields 


o Same descretized form: A . js = b 

o A is a banded matrix, constant 
and same for both <p> and 

o LU decomposition only once 

o Special band solver economizes 
storage and computational effort 

o Judicious implementation of the 
algorithm results in significant 
economy in computer resource 
requirement 
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TURBULENT COMBUSTION MODELING ISSUES 


(FOR GAS TURBINE COMBUSTORS) 


o Most promising method for turbulent reacting flows 

ATTRIBUTES OF DIFFERENT PDF METHODS 


Method 

Joint PDF of <|> 


Joint PDF of 
U and (fc 


Joint PDF of 
U, and co 


Attributes 

Reaction treated 
exactly 


Reaction exact, 
Convection (mean and 
turbulent) exact, 
Variable-density effects 
exact 


Limitations/shortcominas 

Assumes gradient-diffusion, 

Does not give velocityfield 
(requires e.g, k-e) 
Turbulence/chemistry interactions 
not fully simulated 

Needs e equation 
(or equivalent) 


... In addition 

Provides complete closure, 
Treats turbulent streams of 
different scales, 

Can account for effects of 
large scale structures 


PDF CALCULATIONS FOR A RECIRCULATING FLOW 


(Anand et al. 1989) 

x« - (x - x* £ ) /x RE 


*R 



REATTAOWENT LENGTH : 

U rff EXPERIMENTAL (x BC l - 6.75 H 

I X . CALCULATED - 6.65 K 




Backward-facing step 

Pronchick and Kline (1983) 



Storage: 1 .3 Mwords <U>/U ref 

CPU Time: 6 min. Cray YMP 



102<u’2>rtj2 ( 


103< U '2 v ' > /u 3 

ref 
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STOCHASTIC DISSIPATION MODEL 


0 Provides complete closure of the PDF equation (joint velocity-frequency-scalar) 

o More realistic than a mean dissipation model. Dissipation (rather, turbulent 
frequency) is also a random variable and included in the joint PDF. 

o Treats multiple scales in the flow 

o Accounts for internal intermittency 

o Accounts for effects of large scale structures, and influence of origin and history 
of the fluid particles 

! * — — 

J = -CO’ <(D> (So) + C x Q.) dt + <0)> 2 h dt + (0* (2C X <C0> a 2 ) 1/2 dW 

| dU,‘=--i^dt + D i dt + (C 0 kco’) 1/2 dW i 

! P dx i 


SWIRLING FLOWS 


o No theoretical limitations 

o Additional production terms due to non-zero mean swirl velocity 
o Additional terms in calculating the mean pressure (or mean pressure gradients) 

- Boundary layer flows: 

> radial pressure gradient 

> axial pressure gradient also included 

- Elliptic flows 

> additional terms in the Poisson equation for pressure 

o Validation of the stochastic dissipation model and first calculation of swirling 
flows with the joint PDF method (Anand et al. 1 993) 
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/gOINT PDFCC^LC SWI RLiNG FLOWS 


SWIRLING COAXIAL JETS 

(DIFFUSION FLAME SETUP) 
(WPAFB) 


4 UNES - CALCULATIONS 

TL . SYMBOLS -DATA 

; fc ■ ■ JET FLUID 

3 ft ’ • SWIRLER FLUID ! 

v M L ■ - ▲ COFLOW FLUID j 


X/D = 5.29 



V; %.0 02 0.4 0.6 0.8 f 


<u>/<u>„ 



X/D = 5.29 


JOINT PDF CALCULATIONS FOR SWIRLING FLOWS 

COMPARISON WITH REYNOLDS-STRESS MODEL RESULTS AND 
ASSESSMENT OF GRADIENT DIFFUSION MODELING 


x/D = 1.06 



-.40 -.20 0.0 .20 .40 

1 000<u2v>/<U>3oc 


2 _ k 2 3<u > 

<u v> = -C s <V > 

<e> 3r 


C =0.22 


PDF CALCULATIONS 

RS MODEL WITH RSM RESULTS 

RS MODEL WITH PDF RESULTS 
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SPRAY CALCULATIONS 

(Anand 1 990) 


o Advanced spray models (stochastic 
Lagrangian, Monte Carlo) naturally 
compatible with the joint PDF method 

o Assumptions about turbulent kinetic 
energy partition avoided 

o Effects of gas phase turbulence 
structure (velocity cross-correlation) 
included 


./b- 0 S3 1.04 2.87 3.11 6.22 f 2.4 



Computed profiles of normalized 
turbulent kinetic energy of air compared 
against data. 



<U(j>AJ<V C 


«/&- 0.62 1.04 2.07 3. i t 6.22 12.4 



REDUCED KINETICS / MANIFOLD METHODS 

o Low dimensional manifold methods (ILDM, TGLDM) 

- Given detailed kinetes, they provide low-dimensional description 
(e.g., 1-D, 2-D, 3-D) in multidimensional composition/scalar space 

- Use dynmical systems theory to determine the low. dim. manifold 

- Avoid ad hoc assumptions, e.g, partial equilibrium of some of the reactions 
Implications for ignition and lean blow-off 

- Not fuel specific like conventional reduced kinetic schemes 




Perfectly Stirred Reactor (Pope & Maas 1 993) Laminar Premixed Flame (Maas & Pope 1 994 ) 
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PARALLEL PROCESSING 


o Objective: Turnaround time of Iday or less for 3-D combustor calculations 
o Particle partitioning, domain decomposition (multigrid, multi-block) 

o Preliminary results for 2-D flow with particle partitioning (Pope 1994) 

- 16 nodes, 128 MB each, IBM SP1 
- 1 2.8 million particles (800,000 per processor) 

- 50 time steps 

- 44 minutes/processor (45 minutes clock time) 

Extrapolation to 3-D combustor calculations 

- 6.5 hours clock time with 32 processor SP1 


JOINT PDF FOCUS AREAS 


o 3-D Flows, Improved solution algorithms 
o Parallel processing 

o Reduced kinetics / Low Dimensional Manifolds 
o Evaporating / reacting sprays 
o Emphasis on emissions and performance predictions 
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OVERVIEW OF TURBULENCE MODEL DEVELOPMENT AND 

APPLICATIONS AT ROCKETDYNE N95- 27891 


A.H. Hadid, E.D. Lynch, and M.M. Sindir 
Rocketdyne Division 
Rockwell International 
Canoga Park, California 


TURBULENCE MODELING REQUIREMENTS, 
DEVELOPMENT PHILOSOPHY AND APPROACH 

• REQUIREMENTS 

• TURBULENCE MODELING IS A KEY ENABLING TECHNOLOGY 
FOR ALL PROPULSION RELATED CFD ACTIVITIES 

• FACTORS TO CONSIDER INCLUDE ACCURACY, CONSISTENCY, 
COMPUTATIONAL COST, AND EASE OF USE 

• TURBULENCE MODELS THAT CAN NOT BE INCLUDED IN 
PRODUCTION GRADE CFD CODES ARE OF LIMITED VALUE 
TO INDUSTRY 

• PHILOSOPHY 

• BASIC MODEL DEVELOPMENT IS BEST LEFT TO SPECIALIZED 
"CENTERS OF EXCELLENCE" 

• VARIOUS CLASSES OF MODELS NEED TO BE SUPPORTED SINCE 
NO SINGLE UNIVERSAL MODEL IS SHOWN TO EXIST 

• ESTABLISHING THE RANGE OF APPLICABILITY, ACCURACY, AND 
THE COMPUTATIONAL COST OF THE MODELS IS ESSENTIAL 


, TURBULENCE MODELING REQUIREMENTS, 
DEVELOPMENT PHILOSOPHY AND APPROACH (Cont.) 


• APPROACH 


• IDENTIFY KEY "CENTERS OF EXCELLENCE" AND ESTABLISH 
COLLABORATIVE RELATIONSHIP 

• ACQUIRE MODELS AND ASSESS PERFORMANCE FOR THE 
INTENDED CLASS OF APPLICATIONS 

• DELINEATE MODEL DEFICIENCIES AND INITIATE EFFORT TO 
REDUCE THEM 

• DEVELOP MODELS INTO STAND-ALONE MODULES 

• INCLUDE MODULES IN PRODUCTION CODES AND ESTABLISH 
BASELINE FOR APPLICATIONS 
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TWO MAJOR AREAS OF CONCENTRATION 


HIGH SPEED TURBULENCE MODELING (LEAD DR. DOUG LYNCH) 

• FOCUSED ON HIGH SPEED (M>1) PROPULSION (ROCKET 
AND AIRBREATHING) AND AERODYNAMICS 

• EMPHASIS ON 2-EQUATION PHENOMENOLOGICAL MODELS 
WITH NASA ARC AND LARC AS KEY TECHNOLOGY PARTNERS 

• LES WORK IN PLANNING STAGES WITH CTR 


LOW SPEED TURBULENCE MODELING (LEAD DR. ALI HADID) 

• FOCUSED ON LOW SPEED (M<1) AND ROTATING FLOW 
APPLICATIONS 

• EMPHASIS ON REYNOLDS STRESS PHENOMENOLOGICAL 
MODELS IN COLLABORATION WITH UMIST, ICOMP, CTR. AND UAH 

• LES WORK INITIATED WITH CTR 


HIGH SPEED TURBULENCE MODELING 


• EMPHASIS IS ON THE DEVELOPMENT OF ENGINEERING 
TURBULENCE MODELS FOR 

• HIGH SPEED AIRBREATHING PROPULSION SYSTEMS 

• THRUST CHAMBERS 

• VEHICLE AERODYNAMICS 

• APPROACH TAKEN IS BASED ON 2-EQUATION MODELS 

• DIFFERENT CLASSES OF 2-EQUATION MODELS STUDIED 

• k-e 

• k-<o 

• POINTWISE R t 

• COMPRESSIBILITY EFFECTS AND TURBULENCE-CHEMISTRY 
INTERACTIONS MAJOR MODEL UPGRADE THRUSTS 

• COMPRESSIBILITY MODIFICATIONS FROM ARC 

• TURBULENCE-CHEMISTRY INTERACTION MODELS FROM LARC 

• USA AND GASP SERVE AS NUMERICAL PLATFORM 

• GASP - CHIEN. LAM-BREMHORST k-e. k-rn 

• USA - VARIETY OF k-e, k-co 
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COMPRESSIBILITY EFFECTS 

• MIXING LAYER SPREADING REDUCED AT HIGH MACH NUMBERS 

• INCREASE DISSIPATION RATE OF k 

• DEFINE Ck 2 AS A FUNCTION OF TURBULENT MACH NUMBERi^kTyp - 

• ZEMAN MODIFICATION (1990) 

• SARKAR (1990, 1991) AND WILCOX (1991) PROPOSALS 

• MODIFICATIONS OF ZEMAN AND SARKAR NOT RECOMMENDED 
. HEAT TRANSFER OVER PREDICTED NEAR SHOCK WAVES 

• LIMIT TURBULENT LENGTH SCALE L t TO MIN(k^2 ,-H_) 

(VUONG AND COAKLEY, 1987) e 0ji3/4 

. SEPARATION UNDERPREDICTED IN RAPID COMPRESSION OR 
STRAIN REGIONS 

• INCREASE ctjOR a ^UNDER RAPID COMPRESSION (VUONG AND 
COAKLEY) 

■ HEAT TRANSFER OVER PREDICTED FOR VERY COLD WALLS 
Tv/T aw <0.1 (COAKLEY) 

• CEBECI-SMITH ~ 60%, k-(0 - 40%, q -(0 ~ 10%, k-e ~ 30% 


TURBULENCE MODELS ADAPTED TO USA CODE 


ALGEBRAIC 


DAMPING 
WALL LOCAL 


TRANSITION MODEL COMPRESSIBILITY EFFECTS 

BOLMOARY HIGH ORDER MIXING LAYER SEPARATION REATTACHMENT 

CONDITIONS POLYNOMIAL ARNAL SPREADING EXTENT HEAT TRANSFER 


Battwtn-Lomu X 3 Varalons 


X X 


III 

1, Myong-Kasagl 


X X 


k*0 

«*uju2ney2 


1. Sarkar (1991) 4. Yuong & 

2. Zaman (1990) Coaklay 

3. Wilcox (1991) (1987) 


5. Vuong 8 
Coaklay 
(1987) 


2. CMan (1902) X .X k»g 

3. Jonas-Uundar X X 

(1972) t-U 


4. Laundar-Sharma X X 

(1974) 

5. Huaog-Coaklay X X 

(1992) 

6. Sparlala-So-Zhang X X 

(1M3) 

7. laro-Bramhoral X X 

(19*1) 

8. High Ra XX 


ksO 

k>0 

c> 2uik,/y2 
kaO 

ty»0 


1,2,3- «• 

1.2.3. 4. 

2, 3. 4. 

1.2.3. 4. 

1.2.3. 4. 

1,2,3. 4. 


5. 

5. 

5. 

5. 


fca 

1 . High Ra Wilcox 
(1991a) 

2. Low Ra Wilcox 
(1991b) 

q-<p 

Coaklay (1987) 

Ona-Equallon (Goldbarg, - 
Two-TTma Scala1992) 

Ona-Equatton Ry X 

(Goldbarg 1993, 1994) 


k.(l 
k *0 

a ■ T.Zuj/y 2 


k>0 X 

OJy B 0 

kaO 

dy 3 dy ' 


1,2,3. 4. 5. 

1,2,3. 4. 5. 
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M = 3.2 FLAT PLATE FLOW 


• CHIEN k-e MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
COMPRESSIBILITY MODIFICATIONS 

VELOCITY PROFILE TURBULENT VISCOSITY 



REF: G.T. COLEMAN AND J.L. STOLLERY, JFM 56: 741, "HEAT-TRANSFER FROM A 
HYPERSONIC TURBULENT FLOW AT A WEDGE COMPRESSION CORNER" 


MACH 7.05 P LOW OVER AXISYMMETRIC FLARE 


CHIEN k-co MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
COMPRESSIBILITY MODIFICATIONS 


AXISYMMETRIC FLARE 


WALL PRESSURE FOR AXISYMMETRIC FLARE 




REF: Ml KUSSOY AND C.C. HORSTMAN, "DOCUMENTATION OF TWO- AND 
THREE-DIMENSIONAL HYPERSONIC SHOCK-WAVE TURBULENT 
BOUNDARY LAYER INTERACTION FLOW; NASA TM 1-01075. 
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MACH 7.05 FLOW OVER AXISYMMETRIC FLARE 

CHIEN k-a> MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
COMPRESSIBILITY MODIFICATIONS 


WALL HEAT TRANSFER FOR AXISYMMETRIC FLARE 



MACH 8.6 FLOW OVER COLD WALL WEDGE 


M = 8.6 

Tw /Taw - 0.065 



THREE STUDIES 

1. CHIEN k-e MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
CORRECTIONS AND WITH AND WITHOUT MIXING LAYER TREATMENT 

2. HIGH-Re k-co MODEL WITH VARIOUS AIR CHEMISTRY MODELS 

3 BALDWIN-LOMAX TURBULENCE MODEL USING WALL AND LOCAL 
DAMPING 
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MACH 8.6 FLOW OVER COLD WALL WEDGE 

CHIEN k-e MODEL WITH AND WITHOUT MIXING LAYER TREATMENT 


HEAT FLUX CALCULATIONS 



MACH 8.6 FLOW OVER COLD WALL WEDGE 

HIGH-Re k-co MODEL WITH VARIOUS AIR CHEMISTRY MODELS 


HEAT FLUX CALCULATIONS 
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MACH 8.6 FLOW OVER COLD WALL WEDGE 

BALDWIN LOMAX, k-e, k-co MODEL COMPARISONS 


HEAT FLUX CALCULATIONS 



X/L 


LOW SPEED TURBULENCE MODELING 


• EMPHASIS IS ON THE DEVELOPMENT OF ENGINEERING 

TURBULENCE MODELS FOR 

• ROTATING MACHINERY 

• FLOW IN DUCTS AND MANIFOLDS 

• REACTING FLOWS 

• APPROACH TAKEN IS TO 

1 SYSTEMATICALLY ASSESS EXISTING PHENOMENOLOGICAL 
MODELS USING COMMON NAVIER-STOKES SOLVER 

2 IDENTIFY, DEVELOP AND VALIDATE MODEL UPGRADES 
COMMENSURATE WITH OBSERVED FLOWPHYSICS 

3 DEVELOP SELF-CONTAINED TURBULENCE MODEL DECKS 
(MODULES) THAT CAN BE INTEGRATED WITH NAVIER-STOKES 
SOLVERS 

4. PROVIDE GUIDANCE TO EXPERIMENTAL AND THEORETICAL 
RESEARCH IN TURBULENCE MODELING FOR ENGINEERING 
APPLICATIONS 
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TURBULENCE MODELS BEING ASSESSED 


PHENOMENOLOGICAL SINGLE POINT CLOSURE MODELS 



2-EQUATION MODELS ALGEBRAIC REYNOLDS 

K- e (SKEM) STRESS MODELS STRESS MODELS 

(ASM) (RSM) 


NEAR-WALL TREATMENTS INCLUDE (WHERE APPROPRIATE) WALL 
FUNCTIONS, MULTILAYER MODELS, AND LOW-REYNOLDS NUMBER 
APPROXIMATIONS 


TURBULENCE MODEL DECK STRUCTURE AND 
INTEGRATION WITH NAVIER-STOKES SOLVER 
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UAJo 


PROJECT WELL UNDERWAY 


• TEAM 

• MODELS PROVIDED BY UMIST, LERC/ICOMP, ARC/CTR 

• MODULE DEVELOPMENT BY ROCKETDYNE 

• MODULE TESTING BY ROCKETDYNE (REACT, USA) AND 
UAH (MAST) 

• MODEL UPGRADES BY ROCKETDYNE, UMIST, ARC/CTR 

• APPLICATION BY ROCKETDYNE TO TURBOPUMP 
COMPONENT (E.G. IMPELLER) ANALYSIS 

• 2-D MODULES COMPLETED, TESTED, AND RELEASED 

• SINGLE SCALE k-e 

• MULTI SCALE k-e 
■ ASM 

• RSM 


• 3-D MODULE DEVELOPMENT IN PROGRESS 


NONLINEAR ALGEBRAIC-STRESS MODEL 

VORTEX SHEDDING FROM RECTANGULAR CYLINDERS (DURAO, et al) 

PARTICLE STREA KLINES 



0.0 H 1 r — 

i.o o.i o.] 


MEAN AXIAL VELOCITY ALONG CENTERLINE 



. , ' SM 


0-1 0.4 


MEAN KINETIC ENERGY ALONG CENTERLINE 
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ROTATION MODIFIED k-e MODEL 

BACKWARD FACING STEP (DRIVER AND SEEGMILLER) 


STREAMLINE CONTOURS 



p.ir* 


MEAN AXIAL VELOCITY AT X/M=4 


RADIAL TURBULENT INTENSITY (vv/^ref) 




ALGEBRAIC STRESS MODEL 

CONFINED COAXIAL SWIRLING JET FLOW (ROBACK AND JOHNSON) 


GEOMETRY 


STREAMLINE CONTOURS 







REYNOLDS STRESS MODEL (LRR - MODEL) 

BACKWARD FACING STEP (DRIVER AND SEEGMILLER) 

STREAMLINE CONTOURS 



MEAN AXIAL VELOCITY AT X/H=4 AXIAL TURBULENT INTENSITY AT X/H=4 




CONCLUDING REMARKS 

PROGRAMS (BOTH COMMERCIAL AND GOVERNMENT) EMPLOY 
NEW TECHNOLOGY ONLY WHEN IT PROVIDES "ADDED VALUE” 

• REDUCED DEVELOPMENT COST 

• INCREASED RELIABILITY AND PERFORMANCE 

• ENHANCED MANUFACTURABILITY 

THE NEW TECHNOLOGY WE OFFER IS THE COMPUTATIONAL 
ENGINEERING TOOLS FOR PRODUCT DESIGN AND ANALYSIS 

THESE TOOLS ARE THE END PRODUCT FOR ALL ENABLING 
TECHNOLOGY DEVELOPMENT 

• PRE- AND POST PROCESSING 

• ALGORITHMS AND NUMERICAL PLATFORMS 

• PHYSICAL MODELS (E.G. TURBULENCE AND CHEMISTRY) 

FAILURE OF ANY ENABLING TECHNOLOGY JEOPARDIZES THE 
PERFORMANCE (VALUE) OF THE TOOL 

NOW MORE THAN EVER, THERE IS A NEED FOR CLOSER 
COLLABORATION AND COOPERATION BETWEEN GOVERNMENT, 
INDUSTRY, AND RESEARCH INSTITUTIONS TO ENSURE 
MAINTENANCE OF COUNTRY'S TECHNOLOGY BASE 
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RECENT ADVANCES IN PDF MODELING OF TURBULENT REACTING FLOWS 


A.D. Leonard and F. Dai N95“ 27892 

CFD Research Corporation 
Huntsville, Alabama 
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MOTIVATION 

Accurate and Efficient Prediction of Emissions 


1. Accurate Prediction of Emissions From Combustion 
Devices Requires Treatment of Finite-Rate Kinetics 

2. The Effect of Turbulent Fluctuations in Velocity, Energy, 
Composition, etc. on Finite-Rate Chemical Kinetics Must 
be Modeled 
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TURBULENCE/CHEMISTRY INTERACTIONS 
Possible Approaches 


• Neglect Fluctuations 

+ Simple 

- Ignores Effect of Turbulence 

• Eddy Break Up . 

+ Simple 

- Assumes Fast Chemistry 

- Mean Density, Temperature Must Still Be Modeled 

• Prescribed PDF 

+ Efficient 

Limited to Fast Chemistry or Single Step Reaction 

> Composition PDF 

+ Finite-Rate Multi-Step Kinetics 

- Expensive 

- Gradient Diffusion 

• Velocity-Composition PDF 

+ More Accurate 

- More Expensive 


PARTICLE REPRESENTATION 

A Solution Method for a Large Number of 
Independent Variables 


• Computational Requirements Increases Exponentially 
With Dimensions for Finite Difference Methods 


• Computational Requirements Increase Linearly 
With Dimensions for Monte Carlo Methods 
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COMPOSITION PDF SOLUTION 

Stochastic Lagrangian Particle Simulation 

Particle Composition and Position Changed to Model Transport of 
Joint PDF 

• Mean Convection 

- Move Particles Between Cells 


• Chemical Reactions 

- Lookup Table Holds Composition Change 


• Turbulent Diffusion }: 

- Exchange Particles Between Cells LJ- 


• Molecular Mixing 

Particle Interaction Changes Composition 




COUPLING 

PDF Solution is Separate Module 
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CHEMICAL KINETICS 

Reduced Models are Used 

Hydrogen: 2H 2 + 0 2 <=>2H 2 0 

CO: CO + H 2 0 <=>C0 2 + H z 

2H 2 + 0 2 <=>2H 2 0 

Methane: CH 4 + 2H + 2H 2 0 -»CO + 4H 2 

CO + H 2 O~C0 2 + H 2 
2Hj + 0 2 -*2H 2 0 
3H 2 +0 2 »2H 2 0 +2H 


Hydrocarbon:C n +(-^-)0 2 _^ n CO+(n+l)H 2 
CnH^+n H 2 0-*n CO + (2n + im 
C0 + H 2 0 oC0 2 + h 2 
2H 2 + 0 2 <=>2H2 O 


Thermal NO: 


N 2 +0«N0 + N 
N +0 2 «N0 + 0 
N +OH <=> NO + H 


RESULTS TO BE PRESENTED 

• Jet Diffusion Flame (Hydrogen with Helium Dilution) 

• Bluff Body Stabilized Flame (H 2 /CO) 

. Piloted Jet Diffusion Flame (Methane) 

• Generic Gas Turbine Combustor (Propane) 
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HYDROGEN JET DIFFUSION FLAME 


Illustration of Experiment at 
Sandia National Lab 


coubst 

ROW 





Re-10 4 


Fuel 

100% H 2 
80% H 2 , 20% He 
60%.H 2 , 40% He 


60% HYDROGEN FLAME 

Scatter Plots of Mixture Fraction and 
NO Mole Fraction 
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UNDILUTED HYDROGEN FLAME 
Conditional Averged NO Mole Fraction 
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HYDROGEN DIFFUSION FLAME 
Dilution Effects on Emmisions Index 



0 5 10 15 20 25 30 35 40 

Dilution (% He) 
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BLUFF BODY STABILIZED DIFFUSION FLAME 


Illustration of Experiment of Correa and Gulati 


IS cm x IS cm wrindtunncf 
top w*8 



♦ r . 

bottom wad 


BLUFF BODY STABILIZED DIFFUSION FLAME 
Composition PDF Predicts Mean Values as 
well as Velocity-Composition PDF 
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PILOTED JET DIFFUSION FLAME 
Illustration of Experiment of Masri et.al 



Conditions 


Flame 

Fuel jet Velocity 

L 

41 m/s 

B 

48m/s 

M 

55 m/s 


ntmxto 



PILOTED JET DIFFUSION FLAME 
Good Agreement with Experimental Data 


Flam* L 


Flam* B 



Flama M 




x/D = 20 


x/D = 30 
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PI LOTED JET DIFFUSION FLAME 

More Accurate Prediction with Monte Carlo PDF 


No PDF Prescribed PDF 


Monte Carlo PDF 





Flame B at x = 20D 


GENERIC GAS TURBINE COMBUSTOR 
Pratt & Whitney Four-Nozzle Sector 
Combustor Tested at Wright Laboratory 
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MO NTE CARLO PDF COMBUSTOR CALCULATION 
Stochastic Particle Traces 



VERTICAL PLANE THROUGH CENTER OF 
FUEL INJECTOR 


Mean CO Mass Fraction Countours 




RUN TIME AND MEMORY 


3D Combustor Calculation 
(68,000 cells) 


Conventional CFD 


CPU Time 

20 hours 

Memory 

80 MBytes 


Monte Carlo PDF 


CPU Time 

100 hours 

Memory 

120 MBytes 


Paral 

el PDF (Projected) 

CPU Time 

25 hours 

25 hours 

25 hours 

25 hours 

Memory 

30 MBytes 

30 MBytes 

30 MBytes 

30 MBytes 


CPU Time for IBM RS/6000 Model 560 


CONCLUSIONS 


• Monte Carlo PDF Solution Successfully Coupled with 
Existing Finite Volume Code 

- Minor Changes to Finite-Volume Code 

- Can be Coupled with Other Codes 

» PDF Solution Method Applied to Turbulent Reacting 
Flows 

- Good Agreement with Data for 2D Case 

- Demonstration of 3D Elliptic Flow 

• PDF Methods Must be Run on Parallel Machines for 
Practical Use 
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EXPERIENCE WITH TURBULENCE INTERACTION AND TURBULENCE-CHEMISTRY 


MODELS AT FLUENT INC. 

D. Choudhury, S.E. Kim, D.P. Tselepidakis, 
and M. Missaghi 
Fluent Inc. 

Lebanon, New Hampshire 


N95- 27893 


Outline of Talk 


• Part I: Turbulence Modeling 

— Challenges in Turbulence Modeling 

- Desirable Attributes of Turbulence Models 

- Turbulence Models in FLUENT 
— Examples using FLUENT 

• Part II: Combustion Modeling 

— Turbulence- Chemistry Interaction 

- FLUENT Equilibrium Model 

• Concluding Remarks 


Part I: 

Turbulence Modeling and Industrial 

Flows 


• Many industrial flows axe turbulent; certainly in the markets 
that two of our codes, FLUENTand RAMPANT, are focused in. 

• Turbulence augments rates of mass, momentum and heat trans- 
fer, often by orders of magnitude. 

• Most combustion processes involve turbulence and often de- 
pend on it. 

• Choice of turbulence model dictates the accuracy of CFD pre- 
dictions. 

• There is still a large gap between the state-of-the-art and users 1 
expectations and needs. 
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Challenges in Turbulence Modeling 


• Modeling the correlations: puju'- and pu-<£ . 

— Closures based on the “eddy-viscosity” concept (indus- 
try’s most popular choice) 

— Closures based on transport equations (RSM) 

• Modeling an additional transport equation for a scalar quantity 
to fix the state of turbulence. 

- Most popular choice: the kinetic energy dissipation rate, 
e . 

— However, this equation is derived by continuum mechanics- 
based phenomenological considerations and intuition. 

• Modeling of the viscosity-affected, near- wall laminar sublayer. 

— Most popular choice: “Wall-functions” that bridge the 
turbulent field to the solid wall. 

— However, assumptions involved are not always right. 


Desirable Attributes of Turbulence 
Models in Commercial CFD Codes 


• Accuracy and Universality 

— The range of applicability should be as broad as possible. 

— Applicable to complex geometries and unstructured meshes. 

• Economy 

- Mathematically simple. 

— Memory and CPU requirements should be moderate and 
affordable (model formulation and grid distribution re- 
quirements). 

• Robustness 

— Model should be able to solve a wide range of problems 
with little or no convergence problems. 

- Computationally efficient (fast execution speed and uses 
memory sparingly). 
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Turbulence Models in FLUENT 


• k-e model adequate for simple flows with no significant strain 
rates. 

• RNG k-e model for separated flows, flows with large streamline 
curvature, swirling flows, or flows with significant strain rates. 

• RSM recommended for swirling flows or highly anisotropic flows. 


k-e Model: Some Comments 


• Well- tested, used for over 20 years, limitations well understood. 

• It forms a good compromise between universality and economy 
of use for many engineering problems. 

• Subject to the inherent limitations of the Boussinesq’s hy- 
pothesis, i.e., isotropic eddy- viscosity and Newtonian closure 
(gradient-diffusion model). 

• Many assumptions are introduced in deriving the modeled equa- 
tions for the turbulent quantities, particularly the e-equation, 
making their fidelity limited. 

• The constants in the modeled equations are calibrated against 
simple benchmark experiments. 

• As a result, the k-e model performs poorly in flows with cur- 
vature, swirl, rotation, separated flows, low- Reynolds number 
flows, strongly anisotropic flows, etc. 
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Renormalization Group (RNG) 
Based k-e Model 


• Basic theory and derivation are described in Yakhot and Orszag 
(1986). Further details and applications are in Yakhot, Orszag, 
Thangam, Speziale, and Gatski (1992), SpeziaJe and Thangam 
(1992). 

• First introduced in a commercial code, FLUENT, in 1992. 

• The RNG method is essentially a scale-elimination technique 
that can be applicable to the Navier-Stokes and other scalar 
transport equations as well. 

• Removal of successively large scales leads to differential trans- 
port equation models and associated formula for quantities such 
as the turbulent Prandtl/Schmidt number. 

• The basic form of the RNG-based k-e equations remains largely 
the same with the standard k-e model. But, the constants in 
the model equations are derived explicitly from theory. 

• The 6-equation ends up with an additional source term, a 
strain-dependent term. 

• The RNG model can be integrated directly to a solid wall with- 
out using ad hoc damping functions or damping terms used in 
many near- wall models. 


• High- Re form of the turbulence kinetic energy and dissipation 
rate equations derived by RNG procedure are: 


d k dk 

dt + 1 dxi 


= Pk- 




% + U,g- = 1.42 £ft- 1.68 £ + 
at dxi k k dxi \<r e ox t J 


where: 


o’* = <j c — 0.7179 

P k = 2 i*r Si j Si j is the kinetic energy production 


Sij = \ ^ j is the mean rate of strain tensor 


Vt 

» = 


_ q hi 


1 + 0 k 


rj = Sk/e, S = (2 5, 
77o = 4.38, P = 0.015 


ij Sij) 


\ 

2 
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RNG-Based k-e Model (Cont’d) 


• In the low- Re RNG model, a differential relationship exists 
between ^ and u e s (Yakhot and Orszag, 1986). 

• The turbulent Prandtl/Schmidt number is no longer a constant, 
and computed from relationships relating the local value of the 
number to the viscosity ratio (Yakhot and Orszag, 1986). 

• In these relations, as v — * 1, a — ► c*o (the low- Re limit) and 
as v — ► co, g — a -1 — * 0.7179 (the high-Re limit). Here: 

v - Ves/vo, where = i/ 0 + vt 

ct — inverse turbulent Prandti number (a -1 ) 

at 0 = inverse molecular Prandti number (&Q 1 ) 


• In the low- Re regions, <r* and cr € are obtained similarly from 
the Prandti number relationships, with ao = 1.0. 

• The relationships ensure that in the high-Re number part of 
the flow where v » 1: 


= v? = 0.085 — 

and the effective viscosity varies smoothly from the molecular 
viscosity to the turbulent viscosity. 

• The low-Re eddy- viscosity formula does not explicitly involve 
any geometric length scale, i.e., the distance from a solid wall 
used in the damping functions adopted by most low-Re near- 
wall models, which is a very convenient feature for calculations 
for complex three-dimensional geometries. 

• In collaboration with the originators of the RNG model, Drs. 
Yakhot and Orszag, the model has been extended to account 
for the effects of compressibility, swirl, rotation, and premixed 
combustion. 

• The RNG-based k-e model also works well with conventional 
and enhanced (non-equilibrium) wall functions available in Flu- 
ent Inc.codes. 
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The Reynolds-Stress Model in 
FLUENT 


• RSM solves transport equations for the Reynolds stresses: u^uj 
(4 equations in 2D problems, 6 equations in 3D problems). 

• RSM is the level of modeling that has a well established track 
record of out-performing eddy- viscosity models in complex flows. 

• It is computationally more expensive and more inclined to di- 
vergence and stability problems. 

• The simple and most widely tested form of the Launder, Reece 
and Rodi (1975) form is used. 

• The interpolation technique for co-located grids of Rhie and 
Chow (1983) is used. 

• It offers the best choice for highly anisotropic flows. 


Example 1: 

Circle-to-R,ectangle Transition Duct 


• Measured by Davis (1991). 

• Re D =3.9 x 10 5 . 

• Solution Domain. 

— Upstream Inlet Boundary: x/D = —1.0 
— Downstream Exit Boundary: x/D = 8.0 
— A Quadrant of the duct modeled. 



STATION l 
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Axial Velocity (m/s) 


Example 2: 
Cyclone Spar at or 


• Measured by Qing (1983). 

• RSM is used. 

• Cylindrical 55 x 23 x 41 grid. 
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Example 3: 

180° Bend of Square- Cross Section 

• Solution Domain 

- Upstream Boundary: 5.0 D a from the start of the bend 

- Downstream Boundary: 5.0 Da from the end of the bend 
— A symmetric half of the duct modeled. 

• Mesh 

— Orthogonal 101 x 47 x 27 
— Distance from the wall « 0.01£>» 



0.00 0.25 0.50 0.75 1.00 

30 ° 





inside outside 



/ 3 >£° 


Part II: 

Combustion Modeling 


• With environmental awareness, legislations on combustion- gen- 
erated pollutants such as NO x , SO x carbon monoxide, soot, 
unbumt hydrocarbons, etc. have become increasingly tougher. 

• Combustion simulation in industrial applications can help us to 
design combustors with higher efficiencies and lower pollutant 
emissions. 

• The combustion process involves some of the most complex 
phenomena such as chemistry, multiphase flow, turbulence, 
heat transfer and the interaction between these phenomena. 

• Here we focus on gaseous combustion in which the reactants 
may be mixed or non-mixed prior to flowing into the combustor. 


140 




Turbulence-Chemistry Interaction 


• Accurate simulation of turbulent combustion requires a thor- 
ough assessment of the way turbulence and chemistry interact. 
The reaction rate and flame structure primarily depend on this 
interaction. 

• In turbulent flames, chemical rates can be significantly differ- 
ent than those in laminar flames (sometimes several orders of 
magnitude), and the mean chemical rate is not the same as the 
rate calculated based on mean values of the various scalars: 

r( 0 u 0 2 y ...)±r( 8 u 0 2y ...) 

• Turbulent-chemistry interaction is best characterized by the 
Damkohler number which is the ratio of characteristic flow time 
to chemical reaction time: 


Da = 


Z£ — J-Ul 

Tr ~ If V' 


• When Da « 1 chemical reactions axe orders of magnitude 
slower than turbulent mixing and the influence of turbulence 
on reaction can be neglected. 

• When Da >> 1 chemical reactions are very fast and hence 
combustion is controlled by turbulent mixing. 


• At high Da we can exploit the laminar flame concept: turbu- 
lent flame is comprised of an array of laminar flames (flamelets). 
Hence chemical rate expressions can be those obtained in lam- 
inar flames and the effect of turbulence can be characterized 
through the probability density function (pdf). 

• For turbulent diffusion flame, the pdf is usually expressed in 
terms of a scalar which can best characterize mixing, e.g., the 
mixture fraction. Since the rate of reaction is much higher 
than the mixing rate, we can assume that the reaction system 
is at equilibrium. The effect of turbulence is simply felt by the 
fluctuations in the mixture fraction. The mean value of an}' 
scalar in the flame is simply: 

$ = /o 1 mp(.m 

• For turbulent premixed flames the pdf is usually expressed 
in terms of a scalar which can best characterize the reaction 
progress, e.g., normalized temperature: 

f = /q 1 r(c)p(c)dc 
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FLUENT Equilibrium Model 


• For turbulent diffusion flames we use a two-moment beta pdf 
and eq uili bri um data to calculate various thermo-chemical scalars 
in the flame. 

• To obtain equilibrium data we use the popular CHEMKIN li- 
brary of SANDIA, fully interfaced with our codes. CHEMKIN 
contains data on all important gaseous fuels, combustion inter- 
mediates and products as well as their properties. 

• We obtain the mean mixture fraction and its variance from 
their respective conservation equations: 



• To save computational time we calculate the integrals before 
the CFD calculations. 


Concluding Remarks 


• As of now, we provide our users with three turbulence models: 

— the “conventional” k~e model, 

— the ReNormalization Group model, 

— the Reynolds- Stress Model. 

• The Renormalization group fc-e model has broadened the range 
of applicability of two-equation turbulence models. 

• The Reynolds-stress model has proved useful for strongly anisotropic 
flows such as those encountered in cyclones, swirlers and com- 
bustors. 

• Issues remain, such as near- wall closure, with all classes of mod- 
els. 

• Collaborative research with ICOMP will not only serve to fur- 
ther quantify applicability of turbulence models but may bring 
to market new ideas in the field of turbulence modeling for 
industrial flows. 
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EXPERIENCES WITH TWO-EQUATION TURBULENCE MODELS 


Ashok K. Singhal, Yong G. Lai, and Ram K. Awa 
CFD Research Corporation 
Huntsville, Alabama 


N95- 27894 


OUTLINE 


• Introduction to CFDRC 

• Experiences with 2-Equation Models 

- Models Used 

- Numerical Difficulties 

- Validation and Applications 
Strengths & Weaknesses 

• Answers to Three Questions (Posed by Workshop 
Organizing Committee) 

1 . What Are Your Customers Telling You? 

2. What Are You Doing In-House? 

3. How Can NASA-CMOTT Help? 


INTRODUCTION TO CFDRC 


• Young and Energetic (Turbulent) Organization, Dedicated to the 
Continuous Process of Advancement and Effective Transfer of CFD 
Technology 

PROJECTS 

• TWO TYPES OF ACTIVITIES: 
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INTRODUCTION TO CFDRC (Continued) 


• Objective User of Turbulence Models 

(0, 1, and 2 Equation Models, RSM and LES) 

• Humble Developer, e.g. Monte Carlo Joint Scalar PDF 

• Active Participant in Recent Small Eddies of Turbulence, e.g. 

- Stanford Endeavor: "Collaborative Testing of 
Turbulence Models" 1989-1993 

- National Workshops at: NASA MSFC, LeRC/CMOTT, 
etc. 1987-1994 

- ASME/Fluids Engineering Division, Biathlon, Lake 
Tahoe, June 1994 


TWO-EQUATION MODELS USED 

• Standard k-£ Model (Launder & Spalding, 1974) 

• Low-Re k-e Model (Chien, 1982) 

• Extended k-£ Model (Chen & Kim, 1987) 

• Multiscale k-£ Model (Kim & Chen, 1988) 

• RNG-Based k-£ Model (Yakhot et. al. 1993) 

• 2-Layer k-£ Model (Rodi, 1991) 

• k-E 44 " Models 

• k-oo Model (Wilcox, 1991) 

++ Models with Corrections for: Curvature, Rotation, Buoyancy, 
Compressibility, etc. 
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NUMERICAL DIFFICULTIES 


• Positivity of k & e (or co) Is Not Guaranteed in Iterative 
Algorithms 

• Strong Nonlinearity of Source Terms and Coupling Causes 
Numerical Difficulties 

• Inappropriate Specifications of e (or co) at Boundaries or 
in Initial Conditions May Also Cause Divergence 

• Non-orthogonaitiy of Grids Adds to Difficulties 

• Non-smooth Change Over for Two-Layer Model Hinders 
Convergence 


VALIDATIONS PERFORMED 

• Channel and Pipe Flows 

• Backward-Facing Step 

• Turnaround Duct 

• Swirl-Flow Combustor 

• Rotating Disk Cavities 

• Boundary Layers 

• Jets, Wakes, and Mixing Layers 

• Periodic Wakes Behind Bluff Bodies 

Examples of Successes and Failures 

1) Flow Around a Square Cylinder; 2) 180° Square Duct; 3) S-Shaped 
Annular Diffuser; 4) Dump Combustor; 5) Backward Facing Step 
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FLOW AROUND A SQUARE CYLINDER 

Strouhal Number 

Strouhal Number = fH 


Uo 

f = Frequency of Vortex Shedding 

Model/Expt. 

Time 

Period 

Strouhal 

Number 

H = Obstacle Height 

Expt. 

7.25 

0.138 

U 0 = Freestream Velocity 

Standard k-e 

7.1 

0.141 


2-Layer k-e 

7.1 

0.141 

Notes: 

RNG k-e 

7.6 

0.132 


1. Experiments By Durao, Heitor, and Pereira (1988) 

2. Computations with CFD-ACE 

inlet: 78H Upstream; Outlet: 22H Downstream 

Grid: 120x80 

Time Steps: Over 70 Per Time Period 

Rei: Awa, R.K., Singhal, A.K., Lai, Y.Q., "Numerical Simulation Of Periodic 
and 3-Dimensional, Turbulent Flows With CFD-ACE, " ASME Fluid Dynamics 
Conference, Lake Tahoe, NV, June 19-23, 1994. 


FLOW AROUND A SQUARE CYLINDER 
Instantaneous Streamlines 


Mid-Cycle 
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FLOW IN A 180° SQUARE DUCT 



and 3-Dimensional, Turbulent Flows With CFD-ACE,' ASME Fluid Dynamics 
Conference, Lake Tahoe, NV, June 19-23, 1994. 


FLOW IN A 18CP SQUARE DUCT 
Mean Axial Velocity at 0 = 3 
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S-SHAPED ANNULAR DIFFUSER 


* 


* 


O : Data of Steven & Frr(1973) 

: Reynolds Stress Model 

:K-£ Model 




k-e Model and RNG Model Failed to Predict the Correct Location 
of the Maximum Velocity Downstream 

Computations with CFD-ACE; Publication Under Preparation 


Confined Swirling Flow 
for a Dump Combustor 



O □ : Data ofNejad elal.{1989) 
: Reynolds Sfress Model 



* K-e model failed to preserve the vortex core strength 
near center (see x/h=10 & 18 ) 

* Computational results to be presented at 
1994 ASME Winter Annual Meeting (Chicago) 
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BACKWARD-FACING STEP 


Sensitivity to Grid Refinement 





7 a 


L 

3 

inner layer 


1 j * l 

1*9 V» 

,1 

13 

1 i 

Z» 23 


X/H 




Low-Re Model Requires >30 Nodes in the Inter Layer 


Ref. : "Comparative Study of High and Low Reynolds Number 

Versions of k-e Models," RX Awa, C.E. Smith, A.K. 
SInghal, AIAA-90-0246. 


BACKWARD FACING STEP 


2-Layer Model 


80x60 


r oO o 

O o _ 
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> 
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10 15 20 25 
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Computations with CFD-ACE; To Be Published 











EXAMPLE APPLICATIONS 


• Gas Turbine Combustors 

• Liquid Rocket Engines 

• Seals and Bearing Cavities 

• Impellers, Inducers, and Fans 

• IC Engines 

• CFD Reactors 

• External Aerodynamic Flows 

• Plus Many More 


STRENGTHS & WEAKNESSES 


Strengths of 2-Equation Models 

• Numerically Economical 

• Easy to Modify 

• Reasonable Applicability Within Engineering Accuracy 
Weaknesses 

• Use of Wall Functions Requires First Grid Outside the 
Viscous Sublayer. This is Difficult to achieve, a Priori 

• Low-Re Approach Does Not Offer Overall Advantage. 

• Two-Layer Approach Needs More Work (e.g. Smoothing) 

• Reynolds Analogy Inadequate for Heat-Transfer 
Applications. 

• Effect of Surface Roughness on Turbulence. 
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CMOTT/CP QUESTIONS 


1 . What Are Your Customers Telling You? 


2. What Are You Doing In-House? 


3. How Can NASA-CMOTT Help? 


WHAT ARE CUSTOMERS TELLING? 


• PLEASE Don't Confuse Us, 

with Additional Models and False Hopes 

• Conclusions (Confusion) Over Last 15-Years 

Use k~e Model, with Wall Functions 

Wall Functions, Oh No!, Never!! 

Use Low-Re k~e,.: Which One?, How?? (Good Questions) 

k-e Is No Good,; Neglects Non-Isotropicity, etc., etc. 

- Jump on RSM Wagon, Now! 

It Can Take You Anywhere, Eventually!! 

Look How Great is this k~£ ++ 

When and How to Use it? (Good Questions) 

Look How Accurate is this Scheme, No Numerical Diffusion. 
Don't Contaminate the Solutions with Turbulence 
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WHAT IS CFDRC DOING? 


• Using What is Available, in Best Possible Ways 


• Listening to Both Sides (Model Developers and Users) 


• Trying to Resist Peer Pressures 


• Struggling to Find Resources for Mundane Goals Such as 
Developing Guidelines for Correct Use of Turbulence 
Models 


HOW CAN CMOTT HELP? 


• CMOTT Has Been Providing Commendable Service in the 
Very Difficult Subject: Turbulence 

• "Turbulence Subprogram" Should Help Further 

• Additional Effort is Needed In Many Areas, Such As: 

Near Wall Treatment 

Effect of Surface Roughness 
Economical Heat Transfer Model 

Documentation of Experiences in: 

a) Model Robustness(ln Addition to Accuracy) 

b) Model Sensitivity to Grid Distribution and Boundary 
Conditions 

Transition Model (if Possible Suitable for k~e Framework) 
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HOW CAN CMOTT HELP? (Continued) 

• NASA-CMOTT Is One of the Few Groups Sustaining 
Momentum for Turbulence Modeling. 

• It Is In Unique (Privileged) Position for Embracing the 
Challenge of Developing Specific Recommendations 
(Guidelines) For: 

a) Selection of Adequate Models for Different Class 
of Problems 

b) Correct Use of Each Model 

• The Task Is Difficult But Practical 

• Select Fewer Roads, Post Milestones, and Go Further 

• Move An Inch Closer to Users 

CMOTT 

0-1 I I I 1 I I I I I 1 0 

Developers Users 
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PROGRESS IN SIMULATING INDUSTRIAL FLOWS USING TWO-EQUATION MODELS: 

CAN MORE BE ACHIEVED WITH FURTHER RESEARCH? 

N95- 27895 


Vahe Haroutunian 
Fluid Dynamics International, Inc. 
Evanston, Illinois 


BACKGROUND AND OBJECTIVES 


► Two-equation eddy-viscosity models (TEM’s) are the most cost effective 
for the purposes of applied CFD. Give best accuracy vs. cost balance. 

► There is a lot of confusion about true strengths and limitations of TEM's 
especially that of standard k-z model. 

► We have embarked on extensive study of TEM’s over wide range of flows: 

> Identify true strengths and limitations of standard k-z model. 

> Evaluate other TEM’s. 

> Assess emerging models and novel modeling trends. 

> Identify key areas requiring further research. 

► This talk provides brief review of TEM’s from perspective of applied CFD. 

> It provides objective assessment of both well-known and newer models. 

> It compares model predictions from various TEM's with experiments. 

> It identifies sources of modeling error and gives historical perspective of 
their effects on model performance and assessment. 

> It recommends directions for future research on TEM’s. 


REMARK: 


► Many reported poor predictions of TEM’s are primarily due to 
combination of improper choice of near-wall model and over-diffuse 
numerics. 

► TEM performance can be much improved form further research in: 

> Length scale deteimining equation. 

> Advanced (Anisotropic/Nonlinear) Eddy-viscosity models. 
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INTRODUCTION AND BACKGROUND 


O About FDI 

► Over 10 years in business. 

Primary product FIDAP (Fluid Dynamics Analysis Package). 

O About FIDAP 

► First commercial general-purpose finite element CFD program. 

► Models wide range of flows. 

► Over 700 FIDAP licenses worldwide. 

O FIDAP Turbulence Modeling Capabilities 

► Based on two-equation eddy-viscosity models: 

> Standard k-z model (Launder and Spalding). 

> Extended k-z model (Chen and Kim). 

> RNG k-z model (Yakhot, Orszag, Than gam, Gatski and Speziale). 

► Low-Re near-wall modeling based on two-layer approach: 

> Viscous sublayers spanned by single layer of specialized elements. 

> van Driest’s model used in viscous sublayers. 

> Interpolation functions based on universal flow profiles. 

► Latest turbulence modeling enhancements (to appear soon): 

> Anisotropic eddy-viscosity models. 

> Wilcox’s k -© model. 

> Anisotropic version of the standard k-z model. 

O Typical Industrial User 

► Design engineer. 

► Trained in fluid mechanics and heat/mass transfer. 

► F amili ar with range of flows of interest to his/her organization. 

► NOT CFD expert. 

► Little or no background in turbulence modeling. 


O Turbulence Modeling Requirements of Applied CFD Codes 

► Optimal balance of cost and accuracy: 

t> Turbulence modeling overhead of critical concern. 

> Overall accuracy of ± 15% adequate for most cases. 

► Consistent performance over wide range of flows: 

> Heat/mass transfer 

> 2-D and 3-D (Cartesian, axisymmetric) 

> Complex geometries 

> Transient flows 

► Adaptable to a wide range of complex flow physics: 

> Low-Re effects 

> Variable density/compressibility effects 

> Combustion 

> Two-phase 

► Minimum level of user rnput/intervention: 

> No fine tuning model coefficients and/or solution parameters. 

> No physical input other than boundary and/or initial conditions. 

► No geometry dependence: 

> Distance to wall and/or y + dependence. 

► Stable numerical characteristics. 
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TURBULENCE MODELING CONSIDERATIONS 


O Key Modeling Issues 

1- Accurate modeling of mechanisms governing pu'u'j,pufi', pu^ . 

a) Pressure-scrambling 

b) Body forces 

c) Transport effects 

d) Dissipation 

2- Accurate modeling of characteristic turbulent length scales. 

3- Accurate modeling of low-Re near-wall phenomena. 

O Optimal Level of Turbulence Model for Applied CFD 

► Second-Moment Closures (DSMC’s) and (ASMC’s) 

(+) DSMC’s ideally suited to modeling aspects l-a,b,c above, however, 
(-) DSMC’s costly, especially in 3-D in presence of heat/mass transfer. 
(— ) Geometry dependence in current pressure-scrambling models. 

(— ) ASMC*s perform erratically (1-c above not well modeled). 

(-) ASMC’s numerically less stable (stiff equations). 

► Two-Equation Eddy-Viscosity/Diffusivity Models (TEM’s) 

(+) Least costly. 

(+) No geometry dependence (except some low-Re TEM’s). 

(+) Numerically more stable. 

(-) Conventional TEM’s not suitable for modeling effects l-a,l-b,&l-c. 
(+) Room for significant improvement in predicting effects of complex 
strain and anisotropy through the combined use of improved length 
scale equations and advanced eddy-viscosity models. 

(-) Transport effects (1-c), however, cannot be directly predicted. 


LENGTH SCALE DETERMINING EQUATION 


O THE STANDARD 7z-e MODEL 


P 


Dz 

Dt 



£ _ 8 
.C.-G-^Py 


where , 


G = 



and. 


Mr 


d»i , du j 
dxj dx i J 


dXj 


= 0.09, Ci = L44, c 2 = L92, c k = 1.0, a e = 13 

► Remarks on standard k-t model: 

> Use is made of Boussinesq’s “isotropic” viscosity model. 

> Fine scale isotropy is assumed in modeling e equation. 

> Is high-Ke model. Must be used with suitable near-wall sub-model. 

> Many reported poor predictions are due to improper choice of near-wall 
model, mesh density, discretization scheme and boundary conditions. 

> Model predicts much better than commonly believed, if used properly. 

> It does however have its shortcomings in predicting difficult flows 
involving strong anisotropy and/or non-equilibrium effects - it tends 

to be over-difFuse. It predicts flatter flow profiles, shorter recirculating 
zones, and occasionally does not predict subtle separation bubbles. 
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LENGTH SCALE DETERMINING EQUATION 


O THE EXTENDED k-z MODEL OF CHEN AND KIM 

► Employs modified e equation containing extra generation term. 

► Rationale is tliat in addition to turbulence time scale k/z, there is 
further time scale p kIG characterizing response of e to mean strain. 



= 0.09, c x = 1.15, c 2 = 1.9, = 0.25, o k = 0.75, a, = 1.15 

► Remarks on extended h- e model of Chen and Kim: 

> Is high-Re turbulence model. Needs near-wall model. 

t> Gives similar predictions to standard model in equilibrium flows. 

> We find Chen and Kim's (1987) recommended model produce 
predictions that are too under-diffuse in confined flows. 

> We have timed constants c, = 1.35 and c, = 0.05 to improve performance 

> Revised model gives better results for some well-known benchmark 
flows, but improved predictions over standard model are not realized 
consistently. More experience and possibly fine tuning is needed. 


LENGTH SCALE DETERMINING EQUATION 



► Above version is high-Re turbulence model. Needs near-wall model. 

► Most testing of model has been done with simple near-wall model. 

► Our testing of model with more accurate near-wall model indicates 
that RNG model is often under-diffusive in internal flows and can be 
very over-diffusive in some external flows. 

► We have tuned model constants and obtained better overall predictions. 

c M = 0.0865, Cj = 1.45, c 2 = 1.83, o k = 0.8, o t = 1.15, q 0 = 4.618, P = 0.17 

► Revised model gives better results for some well-known benchmark 
flows, but improved predictions over standard model are not realized 
consistently. More experience and possibly fine tuning is needed. 
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LENGTH SCALE DETERMINING EQUATION 


O Additional Remarks on ENG k-z Model: 

► Interesting development though no major breakthrough. 

► Most model constants are predicted from RNG theory. 

► In applying RNG theory it is assumed that turbulence field has very 
wide spectrum and that mertial .avfr-nmge isisotropig. 

► Values of model constants predicted by RNG theory are approximate 
owing to simplifying assumptions made in applying RNG method. 

► Model predictions critically dependent on additional term R. 

► The R term reflects proposed contributions from fine scale ani sotropy . 

► The R term is not derived and modeled using RNG theory. 

► The R term has essential similarities with extra term in £ eq'n 
of extended k-z model of Chen and Kim. 

► Latest model does not predict von Karman constant. 

► The most notable fact about the RNG k-z model of YOTGS is 
that it challenges the notion of fine scale isotropy of turbulence 

> Thus e (and consequently the characteristic turbulent length scale) 
is assumed to be significantly influenced by the fine scale structure. 
These effects are heuris ideally modeled via the time scale ratio x\ . 

> It is interesting to note that the assumption of fine scale anisotropy 
used in modeling R conflicts with notion of a wide and isotropic 
turbulent spectrum used in applying RNG theory to rest of model. 

> It is more likely that the turbulent length scale is influenced strongly 
by large scale anisotropy as characterized by the anisotropy 
tensor 

> Anisotropic eddy-viscosity models can provide estimates of ay which 
can be used to design improved length scale determining eqn’s. 


1 ADVANCED EDDY-VISCOSITY MODELS (Beyond Boussin esq) 

O Anisotropic Eddy Viscosity Models (AEVM’s) 

► There has been renewed emphasis in developing AEVM’s. 

► Lead to better approximations of the normal and shear stresses and 
therefore turbulence anisotropy effects. 

► In addition to more accurate modeling of P u t u j , AEVM’s could 
potentially be used to improve modeling of: 

> Length scale determining eq'n. 

> Generation rate of turbulence energy. 

► Examples of AEVM’s are: 

> Lumely (1970) 

> Speziale (1987) 

> Yo shiz awa (1984), DIA 

> Rubinstein and Baron (1990), RNG 

> Taulbee (1992) and Speziale (1993), derived from DSMC’s 

> Launder (1993) 


► Remarks: 

> Potential of models ha^e been demonstrated using simple tests. 

> Improvements in accuracy often of second-order in magnitude. 

> Not been extensively tested especially for swirling flows. 

> Anisotropic models not yet extended to turbulent heat/mass fluxes. 

> We are presently investigating AEVM’s of Speziale (1987) and 
Launder (1993). 
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O Wall Function Models 


► Produce over-diffuse solutions in off-equilibrium boundary layers. 

> Often fail to predict separation or vortex shedding. 

► Unfortunately still in extensive use in applied CFD codes. 

i 

O Specialized Finite Element Model (FIDAP) 

► Is essentially two-layer model. 

► Avoids fine near-wall mesh via use of one layer of specialized elements. 

► Employs van Driest's low-Re mixing-length model in near-wall layer. 

► Combines low cost of wall function models with accuracy of two-layer 

models. i 

► y + dependence confined to single layer and transparent to user. [ 

Remarks: i 

► Most of historical testing and verification of TEM*s has been done using 
wall functions. The excess diffusion has lead to much confusion in 

assessing TEM's. | 

► Proper assessment of TEM’s requires at least two-layer models. - 

► Wall function approach is simply unacceptable for applied CFD. * 


IMPACT OF DISCRETIZATION ERROR 


O Sources of Discretization Error: 

► Grid refinement (grid convergence). 

► Location of computational boundaries (e.g. f outlet, inlet, entrainment). 

► Choice of discretization scheme in space and time. 


Remarks: 


► Effect of discretization error has received less attention in turbulence 
model development and testing. 

► Most serious source of error results from discretization of advection 
terms (i.e., the upwinding scheme ). 

► Common but dangerous upwinding strategy is used in many CFD codes: 

> Use accurate unbounded scheme in -mean flow equations . 

> Use inaccurate numerically diffuse scheme in turbulence eouations . 

> Overall scheme is stable but often highly diffusive. 

> Most of development and testing of turbulence models has been done 
using above upwinding strategy. 

> In our computations we employ the accurate streamline upwind 
(SU) scheme in both mean and turbulence equations. 

> Even more accurate schemes are available which are based on 
Petrov-Galerkin finite element formulations. 
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I NUMERICAL RESULTS 

O Free Jets 

► Round jet 

► Plane jet 

O Internal Flows with Separation 

► Flow past backward facing step (Kim et. al) 

► Flow past step in channel with diffuser wall (Driver and Seegmiller) 

► Flow in pipe expansion (Szszepura) 

O Transient Flow (V ortex Shedding) 

► Flow past square prism (Lyn) 

O 3-D Flow 

► Flow past passenger car models 


REMARKS: 


► Five sets of model predictions are presented: 

> Standard k-z model 

> Extended k-z model (original) 

> Extended k-z model (revised) 

> RNG k-z model with (original) 

> RNG k-z model with (revised) 


FREE JETS 

The Submerged Plane and Round Jets 



Plane Jet 

Round Jet i 

dSfdx 


wmm 

% error 

Experiment 

=0.105 


«0.095 


Standard k-z model 


-1 

0.112 

18 

Extended k-z model (original) 


-5 

0.10 

5 

Extended k-z model (revised) 


—3 

0.104 

9.5 

RNG k-z model (original) 

0.131 

25 


65 


0.101 

-4 


19 
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TURBULENT FLOW OVER BACKWARD FACING STEP 
Kim et al Test Case: Re = 45000 



x* 


Experiment 

7.0 ±0.5 



6.5 

-7.1 

■Rrtp.nded k-z model (original) 

8.4 

20.0 

Ertended k-z model (revised) 

7.1 

1.4 

RNG k-z model (original) 

7.5 

7.1 

RNG k-z model (revised) 

7.46 

6.6 



TURBULENT FLOW OVER STEP IN CHANNEL WITH 
DIFFUSER WALL 

Driver and Seegmiller Test Case: Re = 3GOOO 



Angle 1 

0 degrees 

6 degrees ! 


% error 


% error 

Experiment 

6.2 


8.1 


Standard k-z model 


-14.5 

EM 

-18.5 

Extended k-z model (original) 

6.6 

6.5 


17.9 


5.76 

-7.1 

WEM 

-8.6 


6.17 

-0.5 

8.33 

2.8 



-1.5 

8.33 

2.8 
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TURBULENT FLOW IN PIPE EXPANSION 


Szczepura Test Case; Re = 890,000 



x* 

% error 


9.51 


Standard k-z model 

9.59 

0.9 

Extended k-z model (original) | 

12.44 

30.8 

EgSESHSEBSgai 

10.6 

11.5 

RNG k-z model (original) 

11.35 

19.5 

RNG k-z model (revised) 

11.39 

20 



TURBULENT FLOW PAST SQUARE PRISM 
Lyn's Test Case; Re = 21400 



Strouhal No. 

D 

Q 

Experiment 

0.132 ±0.035 

= 2.0 

NA 

Standard k-z model 


1.68 

0 

Extended k-z model (original) 

0.131 


-0.07 

Extended k-z model (revised) 

0.135 


0 

RNG k-z model (original) 



-0.07 

RNG k-z model (revised) 

0.133 

BIS 
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1 CONCLUSIONS AND DIRECTIONS FOR FURTHER RESEARCH | 


► For applied CFD, TEM’s strike balance between accuracy and efficiency. 

► The use of inadequate near-wall models and over-diffuse numerical 
schemes obscures true performance characteristics of TEM’s. And this has 
lead to much confusion in evaluation of TEM's. 

► Consequences of using better near- wall model and accurate numerics are: 

> Standard k-z model performs much better than commonly believed. 

> Extended k-z model with original set of model constants produces 
under-diffuse predictions. 

> ENG k-z model with original set of model constants gives predictions 
that can be both under-diffusive or over-diffusive depending on flow. 

> The extended and ENG models with revised set of model constants 
perform better than with original set of model constants. 

► Newer models are quite promising, but do not yet perform consistently 
better than standard k-z model. 

► Significant advances in TEM capabilities may potentially result from 
further research in two key areas: 

> Advanced constitutive-type laws for the Reynolds stresses: 

• AEVM’s appear to be best candidates. 

> Improved length scale determining equation: 

• Better modeling of off-equilibrium effects. 

• Better modeling of large-scale anisotropy effects. 
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TURBULENCE MODELING NEEDS OF COMMERCIAL CFD CODES: COMPLEX FLOWS 


IN THE AEROSPACE AND AUTOMOTIVE INDUSTRIES 


N95- 27896 


Bizhan A. Befrui 
adapco 

Melville, New York 


CONTENT OF PRESENTATION 


. STAR-CD: COMPUTATIONAL FEATURES 

• STAR-CD: TURBULENCE MODELS 

• COMMON FEATURES OF INDUSTRIAL 
COMPLEX FLOWS 

. INDUSTRY-SPECIFIC CFD DEVELOPMENT 
REQUIREMENTS 

• INDUSTRIAL COMPLEX FLOWS: 

APPLICATIONS & EXPERIENCES 

- FLOW IN ROTATING DISC CAVITIES 

- DIFFUSION HOLE FILM COOUNG 
• INTERNAL BLADE COOLING 

- EXTERNAL CAR AERODYNAMICS 

• CONCLUSION: TURBULENCE MODELING 

NEEDS 


STAR-CD: COMPUTATIONAL FEATURES 


• BODY-FITTED NON-ORTHOGONAL 
COORDINATE SYSTEM 

• UNSTRUCTURED COMPUTATIONAL MESH, 
DIFFERENT CELL TOPOLOGIES, IMBEDDED 
MESH REFINEMENT, DISCONTINUOUS 
MESH INTERFACE, MOVING BOUNDARY 
AND INTERNAL INTERFACES 

• PRIMITIVE VARIABLE, SELF-ADAPTIVE 
ELLIPTIC-HYPERBOLIC PRESSURE 
CORRECTION METHOD 

. COLLOCATED-VARIABLE ARRANGEMENT 

• EULER-IMPLICIT TEMPORAL INTEGRATION 

. UD, CD, LUD, SFCD SPATIAL 
DISCRETIZATION, WITH BLENDING 
CAPABILITY 
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STAR-CD: TURBULENCE MODELS 


. TWO-EQUATION MODEL 

- STANDARD k-e WITH CORRECTIONS FOR 
BULK DILATATION AND BUOYANCY 

- HIGH REYNOLDS NO. RNG BASED k-e 
MODEL 

. TWO-ZONE (TWO-LAYER) MODEL 

- HIGH REYNOLDS NO.: k-e VARIANTS 

- LOW REYNOLDS NO.: k -t VARIANTS, 

PRANDTL MIXING 
LENGTH MODEL 


STAR-CD: TURBULENCE MODELS 


. REYNOLDS STRESS TRANSPORT MODEL 

- TRANSPORT EQUATIONS FOR 
CARTESIAN STRESS TENSOR IN NON- 
ORTHOGONAL COORDINATE SYSTEM. 
ON NON-STRUCTURED MESH 

- LAUNDER. RODI, REECE (1975) 
FORMULATION WITH LAUNDER (1989) 
MODEL CONSTANTS 

- GIBSON & LAUNDER (1978) WALL 
REFLECTION MODEL 
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Sep 94 
VIEW 

.000 

.000 

1.000 

ANGLE 

.000 

DISTANCE 

MO 

CENTER 

.076 

.057 

.003 

EHIDDEN PLOT 


Y 

e x 

Driver & Seegmtter Backward Facing Step 
How Domain « -20*H to 3TH 
Mesh * 105 ( Axial) x 45 (Radfal) 



Driver & Seegmtter Backward Facing Step 
Data Inlet B.C. 

Location X/H - 1.5 
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ROTATIONAL SPEED 
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2 S*p94 
PRESSURE 

N/M"2 

LOCAL MX- 329.9 
LOCAL MN- *2937 



0.5 LUO DIFFERENCING 

TURBULENCE MOD& KE - UPPER, RSM - LOWER 


COMMON FEATURES OF INDUSTRIAL 
COMPLEX FLOWS 


. THREE DIMENSIONAL WITH MULTIPLE FLOW 
"COMPLEXITIES" 

- BODY-FORCE FIELDS 

- STREAM SURFACE CURVATURE 

- STRONG PRESSURE GRADIENTS 

- COMPRESSIBILITY EFFECTS 

- LAMINAR-TURBULENT TRANSITION 

- COMBUSTION, SHOCK, MULTIPHASE, NON- 
NEWTONIAN 

. LARGE SCALE DOMAIN AND COMPLEX 
GEOMETRIC CONFIGURATION 

• IRREGULAR. UNSTRUCTURED 
COMPUTATION AL WI SH 

. SPATIAL RESOLUTION DIFFICULT TO ACHIEVE 
ON 0(1 0 s - 10 6 ) MESH CELLS 

. INSUFFICIENT AND UNCERTAIN 
EXPERIMENTAL DATA FOR TURBULENCE 
MODEL VALIDATION/IDENTIFICATION OF 
DEFICIENCIES 
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INDUSTRY-SPECIFIC CFD DEVELOPMENT 
REQUIREMENTS 


• AUTOMOTIVE INDUSTRY 

- EFFICIENT COMPLEX-GEOMETRY, MOVING-BOUNDARY 
CAPABILITIES 

- MEMORY/SOLUTION PERFORMANCE FOR LARGE 
SCALE DOMAIN CFD SIMULATION 

- DIAGNOSTIC/COMPARATTVE EVALUATION 
OBJECTIVES 

- GEOMETRIC FIDELITY AND SPATIAL RESOLUTION ARE 
PRIMARY ACCURACY FACTORS 

• AEROSPACE INDUSTRY 

- REGULAR AND SMALL-SCALE FLOW DOMAIN (BENCH- 
MARK EXPERIMENTAL MODELS) 

- DESIGN/PERFORMANCE OPTIMIZATION OBJECTIVES 

• NUMERICAL AND TURBULENCE MODEL ACCURACY 
IMPORTANT 

• REQUIREMENTS 

• HEAT TRANSFER 

• LOW REYNOLDS NO. FLOW 

• BODY FORCE FIELDS 
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EXTERIOR TUBE 
(RIGHT HALF REMOVED) 


ORIGINAL PAGE IS 
OF POOR QUALfTY 


SYMMETRY PLANE 


OUTLET 



MOVUG WALL 
SIMULATION OF GROUND 
MOTION RELATIVE TO 
VEHICLE (40 kph) 


simOu^onof rotating tires 

(371 RPM) 


INLET 

SIMULATION OF VEHICLE 
TRAVELING AT 40 kph 
(T «30C, P « 100 kPa) 


FIGURE 1 : EXTERIOR BOUNDARY CONDITIONS FOR W202 40 kph ANALYSIS 
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APPLICATIONS & EXPERIENCES 


APPLICATION 

(DATA) 

FLOW 

COMPLEXITY 

TURBULENCE 

MODEL 

FINDINGS 

T.M. NEEDS 

ROTATING 
DISC CAVITY 1 

• FORCE FIELD 

• WALL 
EFFECT 

• k-e 

• 2 LAYER k-* 

• EKMAN LAYER 
RESOLVED 

• FAIR 
PRESSURE 
DROP 

• EXCESSIVE 
E.V. 

• RSTM + 
SUITABLE 
2 LAYER 

• LOW Re 
RSTM 

DIFFUSION 
HOLE FILM 
COOLING 2 

• JET-CROSS 
FLOW 

• WALL 
ANISOTROPY 

• k-e 

• RNG, k-e 

• 2 LAYER k-( 

• JET 

SEPARATION 
SENSITIVE TO 
MESH 

TOPOLOGY/ 

RESOLUTION 

• POOR 
SPANWISE 
SPREAD 

• RSTM + 
SUITABLE 
2 LAYER 

• LOW Re 
RSTM 


’GRABER et al (1987) 

2 GOLDSTEIN et al (1968), LIGRANI et al (1992) 


COMPRESSOR DRUM TEST RIG STAR-CD CONJUGATE HEAT TRANSFER MODEL 



Axis of Rotation - 
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CFD Dlcrete Hole FUm Coorsng Verification Study 

Simulation ol oxpotlmant ol Goldstein, at. «l. [19681 ; Blot-ng ratio M -55 

Temperature contours on spanwtse planes ; 2-Iayer model. 


184 









APPLICATIONS & EXPERIENCES (cont'd) 


APPLICATION 

(DATA) 

FLOW 

COMPLEXITY 

TURBULENCE 

MODEL 

FINDINGS 

T.M. NEEDS 

INTERNAL 

BLADE 

COOLING 3 

• FORCE FIELD 

• B.L. 

DISRUPTION 

• k-e 

• DEPENDENCE 
ON MESH 
RESOLUTION 

• GOOD AP, h 


EXTERNAL 
CAR AERO- 
DYNAMICS 4 

• B.L. 

STRUCTURE 

INTERACTION 

• COMPLEX 
WAKE 

• k-e 

• RNG k-e 

• 2 LAYER k-/ 

• DEPENDENCE 
ON MESH 
RESOLUTION 

• GOOD C D 

• POOR LIFT 

• RSTM 

• LOW Re 
RSTM 


3 GE AIRCRAFT ENGINES [ABUAF & KERCHER (1991)] 

4 10 FORD 1/4 SCALE MODELS IN WIND TUNNEL TEST [WILLIAMS et al (1994)] 















29 Aug 94 

K T, COEFFICIENTS 


W/SaMETHR-K 


LOCAL MX- .1092E+05 
LOCAL MN- jOOOOE+OO 













experiment results 

COMPARISON OF EXPERIMENTAL AND COMPtmONAL LIFT COEFFICIENTS 
K EPSILON TURBULENCE MODEL * ~ INHIAL RESULTS *" 
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CONCLUSIONS: TURBULENCE MODELING 
I IMMEDIATE NEEDS 


. NEAR-WALL TURBULENCE 

- ECONOMICAL, ROBUST LOW REYNOLDS 
NUMBER 2 EQ. EVM's AND RSTM 

- A GENERAL AND VERSATILE NEAR-WALL 
TREATMENT FOR RSTM 

• RSTM MODEL 

- ALTERNATIVE CLOSURE OF THE WALL 
REFLECTION COMPONENT, WITHO UT NEED OF 
WALL TOPOGRAPHY PARAMETERS 

. EDDY-VISCOSITY MODELS 

- EXTENSION OF THE NON-LINEAR k-e TO 
INCORPORATE FORCE-FIELD EFFECTS 

• BENCHMARKING 

- A RELIABLE DATABASE OF BENCHMARK SET OF 
REPRESENTATIVE COMPLEX FLOWS 

- BENCHMARK PERFORMANCE CLASSIFICATION 
OF VARIOUS EVM’s (k-e, k-®, RNG AND NON- 
LINEAR k-e, MULTISCALE EVM's) AND RSTM 
CLOSURE VARIANTS 


CONCLUSIONS: TURBULENCE MODELING 
PROGRAM NEEDS 


• A LARGER VIEW OF THE RSTM DEVELOPMENT 
TOWARDS IMPLEMENTATION IN GENERAL 
COORDINATE, COMPLEX GEOMETRY DOMAIN, 
UNSTRUCTURED CFD METHOD 

. A BROADER APPLICATION OF DNS TO 
COMPLEX FLOWS TO ASSIST TURBULENCE 
MODEL DEVELOPMENT/OPTIMIZATION 

• WELL-POSED EXPERIMENTAL DATA, 
OBTAINED IN THE ORIGINAL OR REDUCED 
SCALE MODEL OF THE INDUSTRIAL 
COMPONENT FOR CFD VALIDATION 

. COLLABORATIVE INDUSTRY-CFD 
RESEARCH/DEVELOPMENT PROGRAMS FOR 
EXPERIMENTATION - CFD VALIDATION 
(CALIBRATION) FOR SPECIFIC INDUSTRIAL 
APPLICATIONS 
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TURBULENCE REQUIREMENTS OF A COMMERCIAL CFD CODE 


J.P. Van Doormaal, C.M. Mueller, and MJ. Raw N95- 27897 

Advanced Scientific Computing Ltd. 

Waterloo, Ontario, Canada 


Outline 


• Profiles 

- ASC 

- Application 

- Client 

• Needs 

- Clients’ 

- ASC’s 

• ASC Directions 

- Research 

- Development 

- Products 

• How Can CMOTT Help? 
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Profile of ASC 


• Established in 1985 

• Components of business 

- development 

- applications 

- licensing and service 

• Geographic markets 

- North America 

- Europe 

- Pacific rim countries 


Application Profile 

• Rotating machinery components 

- hydraulic turbines 

- pump 

- compressors 

- turbines 

- stators 

- wicket gates 

- scrolls 

- volutes 

- inlets and diffusers 

- seals 

- stage 

- rotor stator 
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Application Profile cont’d 


• Combustion 

- gas turbine combustor 

- coal fired boilers 

- gasification 

- fire suppression 

- emissions reduction 

- safety 

• High speed external - ballistics 

- explosively formed projectiles 

- finned projectiles 

- sabot discard 

• Heat transfer 

- turbine cooling 

- nuclear reactors 

- heat exchangers 

- electronics system cooling 

• Typical uncertainties 

- geometry 

- initial and boundary conditions 

- transient effects 

- transition 

- limitations of physical models 

- numerical error 
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Client Profile 


• Companies or divisions 

- industrial/manufacturing/research 

- 10 - 200 employees 

- limited or no access to high 
performance computing 

• Users 

- design and/or analysis 

- <3 people 

- network of engineering workstations 

- turnaround time in less than a day 
for analysis, hours for design 


Clients’ Needs 


Needs are most readily identified through typical 
questions from clients. 

• General 

- I am using k-s or two-layer or k-to, or RNG .... what does it 
mean to my calculation? Tell me in words what the deficiencies 
of the model means for my application? 

- What is the relative price/performance of the various turbulence 
models? 

- Has the model I am using been validated for type of flows I am 
trying to model? If so, when, where, how ... ? 

- How well does the model handle the interaction between 

turbulence and rotation, curvature, adverse pressure gradients, 
separation, swirl, bouyancy, extinction, droplets and particles, 
anisotropies ...? 

- How can I use Navier-Stokes solvers for design? Can I tune the 
turbulence model to suit my needs? If so, what are the 
appropriate settings for my application? 
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Clients’ Needs cont’d 


. Grid 

- I don’t have access to high performance computing, I 
don’t have any more time, I have a coarse non- 
orthogonal mesh, is my CFD result useful? 

- I have just made my grid finer, why should I have to 
worry about whether y+ is in a given range? 

• High speed flows 

- I am solving a flow with many speed regimes including 
low speed separations and shocks, why do turbulence 
levels become unphysical as the grid is refined through 
shocks? 

- How should experimental data be compared to results 
from time or Favre averaged calculations? 


Clients’ Needs cont’d 


• Combustion 

- Which of the many different combustion models in 
combination with which turbulence model works best for 
my application? 

- How appropriate is the single scale implicit in the 
turbulence model for the combustion model? 

- How can the Bousinesq assumption be valid in the 
presence of counter-gradient diffusion? 

- How important are turbulent fluctuations to my problem? 

- If I had all the mean flow and fluctuating components of 
the the turbulent flow, how can the effects of stretch and 
curvature on the instantaneous flame front be modelled. 

- Can extinction due to vortex stretching be modelled? 

- What is the influence of the flame front on the turbulence? 
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Clients’ Needs cont’d 


• Calculated pdf models 

- If I use a more detailed chemistry model - like a pdf 
transport model - how much improvement can I expect in 
the results for my application? How can I measure that? 

- Is it the case that the results for my application will not be 
sensitive to the shape of the pdf? If not, then why should I 
incur the costs associated with a pdf transport equation. 

- I am solving a pdf transport equation, how much are the 
results dominated by the limitations of modelling of the 
diffusion transport term? 


Clients’ Needs cont’d 


• Flamelet models 

- I am using a flamelet model in modelling my gas turbine 
combustor, but in some regions of the combustor the 
model is not strictly appropriate - can any of the results 
be used? If so, how much? 

- In some models like the flamelet model, it is assumed 
that the turbulent time scale is inversely proportional to 
the velocity gradient of a “laminar” model flame. What is 
the validity of this assumption? 

- How sensitive are my results to the assumption of 
statistical independence of the quantities in a joint pdf? 
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ASC’s Directions 


• Develop in-house model expertise 

- two-layer model 

- alternative two-equation models 

- second moment closure models 

- expanded EBU models 

- flamelet model 

• Develop in-house expertise applying models 

- turbomachinery 

- combustion 

- heat transfer 

• Promote high performance computing 

- parallel computing 


How Can CMOTT Help? 


• Model improvements to 
address between 
turbulence and 

- rotation 

- curvature 

- adverse pressure 
gradients 

- separation 

- swirl 

- bouyancy 

- droplets and particles 

- anisotropies ... 


as well issues related to 

- extinction 

- trace species 

- vortex stretching 

- flame fronts 

- time and length 
scales 

• Great, but is this what 
users really want? 
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How Can CMOTT Help? cont’d 


• Curator of information on existing 
models 

- define 

- validate 

- process 

- educate 

as an independent agency 


How Can CMOTT Help? cont’d 


Define models 

- unified conceptual framework 

- establish baseline for various 
models 

- set context for model improvements 

- for each model 

> document derivation 

> identify assumptions 

> clearly state implications of assumptions 

> separate physics from numerics 
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How Can CMOTT Help? cont’d 


Validate models 

- fundamental flows 

> validate assumptions 

- benchmark problems 

> select real engineering problems 
relevant to identified applications (in 
propulsion) 

> review selection of benchmark on 
regular basis 

- experimental data 

> collect and review existing data 

> define new experiments 

> review quality of resulting data for 
validation of models 


How Can CMOTT Help? cont’d 


Process data 

- collect 

- distil 

- review 

- interpret 

- describe 

- compile 
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How Can CMOTT Help? cont’d 


Educate 

- document 

- publish 

- workshops 

- seminars 

- short courses 

- market 


Summary 


Provide information so users, for their 
applications can: 

• make an educated choice of model 

. understand how to appropriately use 
existing models 

• move forward with existing models and 
technology 

• understand implications of 
improvements to existing models 
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SECOND-ORDER CLOSURES FOR COMPRESSIBLE TURBULENCE 


J.L. Lumley, S. Savarese, and C.C. Volte N95- 27898 

Mechanical and Aerospace Engineering Department 
Cornell University 
Ithaca, New York 


OUTLINE 


I. Project Description 


II. Turbulence Modeling 




III. Computational Engine / Results 


FUTURE WORK 


I. PROJECT DESCRIPTION 


1. Flows of Interest 


2. Motivation 


3. Method 













• Physics 


Boundary Layer Separation & Wall Heat Transfer 
— Spreading rate 

• Modeling 

Account for Compressibility Effects 
on Turbulence 

• Numerics 

Compare 1-point Closures on Identical Solver 
1.3. METHOD 

• 1-Point Closures: from EVM to Second-Order Closures 

• Dynamical Compressibility Effects 

• 3D / Finite Volume Approach 
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II. TURBULENCE MODELING 



1. ("insure Levels 

2. Compressibility Effects 

3 shock Wave I nteractions 
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II- 1. Closure Levels 


1. EVM Mixing-Length 
(Baldwin-Lomax) 


2. EVM Multi-Equation 
(k-e-S) 


3. Second-Order Closure 
(Shih and Lumley) 


II. 2. Compressibility Effects 
1. New Physics & Averaging 


2. Models 



i I MIL 


11.2.1. New Physics 
(Turbulent Kinetic Energy Sink) 

- < TijUiJ >— rid — £ d ~ £ s 

• £d — (MB + jM) < d 2 > 


• n«i =< pd > 


II. 2. 3. Turbulence Modeling 
(Zeman, Sarkar et al. , Yoshizawa) 


• dilatation dissipation: 

£d = (mb + ^m) < d 2 > 

- Sarkar et al. (asymptotic analysis) 

- Zeman 
(Shockiet model) 
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• pressure-dilatation correlation- 


=< pd > 

— Zeman (acoustic model): 

- Sarkar et al. 

(DNS & asymptotic analysis) 
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TT.3- Sho<~k Wave Interactions 


1. Experimental Observations 

2. Physics 

3. Modeling 


11.3.1. Experimental Results 

• Oscillation increases with Shock Strength 
(Dolling) 

• Oscillation increases with Separation Region 

• Normal Stresses Preferentially Amplified 
(D^lery et al.) 
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II. 3. 2. Physics 


Oscillation Caused by (?): 


• "Breathing" of Separation Region 


• Vortex Bursting 
in Incoming Boundary Layer 
(Dolling) 


II. 3. 3. Shock Oscillation Modeling 


• Parametrized Source Terms 

in Normal RS Evolution Equation 
(gradient activated) 

• Separation region Extend 
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HI. COMPUTATIONAL ENGINE 


1. Numerical Method 

2. Turbulence Models 

3. Validation Procedure / Results 


111.1. Numerical Method 


1 Initial Code: flo103 ~| 

Current Code: cvste 

Future 

HIM 1 

(D.Cauqhevl 

1. Geometrv 

Princeton) 

1. Geometry 

3D 

1. Geometry 

0- R-meshes 

2. Turbulence Models 

C-mesh 

(EAGLEView MSU) 

soc 

2D 

2. PDE Solver 


2. PDE Solver 

variable number of PDFs 


spatial discretization: FV 

consistent gradient comp. 


time integration: RK 

3. Convergence Acceleration 


3. Convergence 

Enhanced multigrid sequencing 


Acceleration: 



variable time step 

4. I/O 


residua! smoothing 

Restart option 


artificial dissipation 

Post-processing (DX.Tecpiot,...) 


multigrid preconditioning 

convergence histories 


4. I/O 

5. Turbulence Models 


PLOT3D format 

k-epsilon (-S) 



6. Software Engineering 


5. Turbulence Models 

Dynamical mem. allocation (C) 


Baldwin-Lomax 

Vectorized data structure 



Unix Integration 
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III.2. Turbulence Models: 

Incompressible / Compressible: an additive approach 


• Baldwin-Lomax 

• k-Epsilon / k-Epsilon-S: B.C’s 

• Second-Order Closures 
Boundary Conditions: Wall-Functions 


III-3. Validation Procedure / Results 

• Calibration against simple well-documented 
flows (flat plate, jet ) 

• Results and Comparison of models 
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— 2D => 3D 

— More Complex Wall Functions 

— Realizability Conditions (SOC) 


• Modeling 

— Refinement of Existing Models , < pd>) 
- Shock Oscillation Model 
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ORIGINAL page is 

OF POOR QUALITY 
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Homogeneous Turbulence (R) / k-eps-S / Mach=0.045 Re=24357 


Downstream Evolution 



Homogeneous Turbulence (R) / RSC / Mach-0.5 


Downstream Evolution 




ORIGINAL PAGE IS 
OF POOR QUALITY 
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N95- 27899 


MODELING OF TURBULENT CHEMICAL REACTION 


J.-Y. Chen 

Department of Mechanical Engineering 
University of California, Berkeley 
Berkeley, California 



Modeling Turbulent Reacting Flows 



\ 

Model for 
Chemical Kinetics 
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Regimes of Turbulent Combustion 



Rt 


Turbulent Reactive Flows edited by P A Libby and F A Williams (1994) 


Regimes of Premixed Turbulent 
Combustion 



Turbulent Reactive Flows edited by P A Libby and F A. Williams (1994) 
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Regimes of Non-Premixed Turbulent 
Combustion 



Turbuient Reactive Flows edited by P A. Libby and F.A. Williams ( 1 994) 


Chemical Closure Models 

(1) Laminar Chemistry 

< pWj >=pw i (Y i ,T) 


(2) Fast Chemistry 


1_„ 3 2 Y e (f) 

<PW ' >= - 2 PZ '^^ 


(3) Flamelet model 

< P w i >= JJ P w i (Tl,Xf ) p f. xf (Tl.Ef )dTide f 


(4) Assumed PDF: 

<P w i >= J.Jpw i (<)) i ) P $ d(t) 1 d(l) 2 ..d(f) n 

Assumed the shape of P<h 

(5) Scalar PDF method: 

Solve for Pp directly. 

(6) Conditional Moment Closure (CMC) 

<pwj >=|<pw ; Irj > P f (ri)dr| 
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Flameiet library with one side being burned premixed flamed- 1 .4 



Flamelet Model: 69%H2 +31%CH4 



Mixture Frodlon 


Vraftos, « al. "Nitric Oxide Formation and Differential Diffusion in a Turbulent Methane- 
Hydrogen Diffusion Flame,” 24th Sympostum(!ntcmational) on Combustion/Thc Combustion 
Institute, 1992/pp 377-384 
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Flamelet M odel: 69%H2±3j_%CH4 
Turbu lent Jet Flame. Rev =10.000 



Hydrogen D.ffiision Flame.” 
Institute, 1992 /pp 377-384 


J> Symposiumfljilcmalional) on Combusbon/Thc ConjHisuon 


Flamel et Model: 69%H2+ 3 1 %CH4 
Turbulent Jet Flame ftev.=H) Qon 



vranos. ei al 


„ j Nltnc 0:udc Fofrnalion and Differential Diffusion in a Turbid 

.ns^.^pTj"^" ** Symp0! ' Um(lM ~') « Combustion/ 


Methane- 

Combustion 
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3 extinction properties 


Advanced F^rnelet Approach 


Local models — j 

[ 

Complex kinetics 

Ignition or flame 
stabilization submodels 


Multicomponent 

transport 



flamelet model 
dynamic equations 
specks conservation 
equations 


ill 

S 


Reaction rates, 
extinction conditions 


■ 

_ 


Extinction 

conditions 


l 


Local consumption 

rates 


Balance equation for the flame 
surface density 


Effective strain 
rate model 



■ — 

Flame interaction 

Flame/wall 

model 

interaction 

i 


Vracttons 


Non-flamelet 

combustion 

submodel 




Consumption 

rates 

Heat release 
■ate 


Turbulent flow description 
dynamic equations 
specks conservation equations 
closure rules 


Global model 


r’nnfHtional Momen t Closure (CMC_j 

Definition: 


< Y, lr| >s< Y (x,t)lf(x,t) = 'n> 


Equation: 


<plT , > i5i!2H + <puir,> v<Y i in> + 


at 


V -{< pu'y' ,T l > P f (tj)) 

p f en) 


a 2 <Y l iTi> 

=< pw,iri > + < pDjVf ■ Vflr) > 


Modeling: 

< Wjlri >« Wj(< Tlq >,< Yjlri >.-•) 

< pD.Vf • VflT| >=< pD.Vf -Vf >=^PZt 

< pulr| >= pu 

< pu'y , ' T l >= 0 

< pit) >= p(< Yjltj >,< Tlrj >) 
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.Local strain rate, . 
Equivalence ratio... 







1 5 6 7 

log(Cr,/D,) (l/ s ) 
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Conditional Moment Closure (CMC) 

Applications: 

• Incorporated into existing moment 
closure CFD codes for complex 
geometry flows 

• Realistic Chemistry - Detailed or 
reduced 

Research issues: 

• Modeling of conditional statistics 

• Preferential diffusion 

• Parallel computing algorithms 


Probability Density Function (PDF) 

Applications: 

. NO x from methane jet flames with 
reduced chemistry 

• Sooting flames 

• 2-D flows 

Research Topics: 

• Mixing model 

• Extension to droplet spray & particle 
laden flows 

• Preferential diffusion 

• Efficient stochastic algorithm 

• Construction of chemical tables 

• Parallel computing - 3D Flows or 2D 
flows with complex chemistry 
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Mixing Models for PDF Methods 


. Modified Curl's Model (stochastic) 

v d 2 J/e i5=vjAp_(y,0 


(I, 

V V 


/L(y\r)P-(vy V)#(V'. V"l'i')- V V,,) ^ W 
d) t 


- IEM (Interaction-by-Exchange-with-the-Mean) Model 

(deterministic) 


Mixing Frequency: <* mix = ft mix 


c » a 


' 2T ma^a 


(V-«^CV.o| 


Modified Curt s Model 




PaSR: H2/NOx Detailed Chemistry 4> 1 t lm£ 


1.0x10' 


0.8X10 i 


£ 0.6X10 


5 0.4x10* H 


O 

z 


0.2x10" 


• modified Curt's model 
a lem mixing model 


0.1 02 
Unmixed ness 
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Comparison of Predicted and Measured 
H20 Mass Fractions 
Turbulent Nonpremixed Jet Flames 


Fuel: (CO/H2/N2: 0.30/0.10/0.6 ) 



Experim ental Evidence of Preferential Diffusion 
in Turbulent Jet Flames 
(Fuel: 36%H 2 +64% C0 2 ) 



Differential Molecular Diffusion In Reacting and Nonreacting Turbulent Jets of H2/C02 
mixing with Air. ' LL. Smith Ph D Thesis. University of California at Berkeley (1994) 
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INTRODUCTION TO TURBULENCE SUBPROGRAM 
T.-H. Shih and J. Zhu 

Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 
NASA Lewis Research Center 
Cleveland, Ohio 


OBJECTIVES 


• A means for CMOTT to interact with industry 

• A vehicle for technology transfer to industry 


CONCEPT OF TURBULENCE MODULE 
Exact CFD equations: 


_ d dUi dUj 2 dU k 


Dt 


dP 


Reynolds stresses will be recasted as: 


■dUi dUj 2 dU k 


—pUiUj = u T ( — 1 j- zsLl w uu k f . r qu ATT o 

= ^ 


*i)j 




CFD equations become: 


D m - 4 ^ + + U - 5^)1 + g f 

• The task of turbulence module: Provide ^ and t 
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• Turbulence Module: 

0 Input: Ui, P and m • • ■ from the mean flow solver 
0 Output: 


r fc2 f^ = 

VT = C»— [ Dt 

.dUi^dUj 

T 'i = -f> u ' u 3 -^ T{ d^ + dx t 


De 

Dt 


— 1 


2 31/t , , 

iS * 1 ’ 1 


0 Models for p muj 

- One- and two-equation eddy viscosity models 

- Reynolds stress algebraic equation models 

- Reynolds stress transport equation models 
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Module with CMOTT research code (incompressible) 

• CFD equations in CMOTT research code: 


DpUi 

Dt 


d u ..dUi - , 


dUj 


dP 

dxi 


• Turbulence module: provide p. T and T tJ 
0 Built-in models without wall function: 

- Launder-Sharma and Chien k - e models 

- CMOTT k-e model 

0 Built-in models with wall function: 

- k-ui model, standard k-e model 

- CMOTT k-e model 

- CMOTT Reynolds stress algebraic equation model 


Module with NPARC code 

• CFD equations in NPARC code: 


DpUi 

Dt 


dr dUi dUj 


dxj 


dx. 


dx, 


2dU k dP_ 

3 dx k lj )i dx { 


• Turbulence module (present time): provide isotropic ht 
0 Build-in models without wall function: 

- Baldwin-Lomax model and Chien k-e model 

- CMOTT k-e model 
0 Further development: 

- Models with wall function 

- Reynolds stress algebraic equation models 

- Reynolds stress transport equation models 
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Joint Program with Industry 
on Turbulence Module 

• For those who want to use the available modules: 

0 Need interface program for particular industry codes 
- Grid informations, Boundary treatment, ... 

• For those who want a module for their own codes: 

0 Need modules exclusively for particular industry codes 

• Maintain and update the turbulence modules along with 
model development. 
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DESCRIPTION OF TURBULENCE SUB-PROGRAM 


J. Zhu 

Institute for Computational Mechanics in Propulsion 
NASA Lewis Research Center 
Cleveland, Ohio 


General Transport Equations 

Dit ) = rJ-' Sl) , 

• Non-dimensional form ( M , <* ^/ Ret w / fle ) 

• Conservative form 

• Cartesian velocity components 

1. Easy to transform (chain rule) 

2. No curvature terms 
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Discretization 

• Finite-volume method 

• Source term 

S^ = Si + S 2 <t>, Si>OandS 2 <0 

• Transient term 

1. lst-order fully implicit scheme 

2. 2 nd-order three-level fully implicit scheme 

• Diffusion term 

Standard central differencing scheme 


• Convection term: HLPA scheme 

(Hybrid Linear/Parabolic Approximation) 

^ _ 4>w - <Pww 

(j) w — (fiw + 7 (0C — — (fiww 


7 = 



if \4> w - 0.5] < 0.5 
otherwise 


- second-order accurate 

- Bounded (non-oscilletory) 

- Diagonally dominant matrix 


•N 



s 

•S 


•SS 
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Solution Procedure 


• Non-delta form 

Positiveness (</> > 0 but A0 may < 0) 

Simple linearization 

• Algebraic equations 

A c 4>c = A w 4> w + A e 4>e + A s<t>s + a n4>n + s 
A's, S> 0 

• Decoupled solution 

• Alternating direction TDMA solver 


Boundary Conditions 

• Inflow: <p specified 

• Outflow: Fully-developed condition 

• Symmetry: dfy/dn = 0 

• Wall: 

1. Low-Reynolds number turbulence models 

2. Standard wall-function approach 
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Sub-Programs 

• NPARC2D version 

Plane or axisymmetric, without swirling 

Compressible 

Non-vectorized 

• FAST2D version 

Plane or axisymmetric, with or without swirling 

Incompressible 

Vectorized 


IMPARC2D Version 

• Grid arrangement 
Control volume centers 
Boundary nodes 
Embedded bodies 




/// 

/// 

IZZZZZZZZZZZZZ 


UT77T7TT77\ 


J-Patches 



J 



K-Patches 
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• Input from the main code 

1. Geometric quantities: x, y, £®, Vx, ly, J 

2. Flow variables: /x, J _1 p, J^pU, J 1 pV, J 1 E 

3. Patch control: 5x2 parameters 

4. Boundary conditions: 7x2 parameters 

• Output 

1. To the main code: yt 

2. For post-processing: K , e, y + , 2 /n, fn 


FAST2D Version 


• Grid arrangement 
CV centers 

B1 

Boundary nodes 
Embedded bodies 
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• Vectorization 
Single-index: 
ii=i+(j-l)ni 
<t>0 J)=^(ii) 

0(iJ-l)=0(ii-ni) 
Control parameter: 



KBLK = I 1 for computational nodes 
\ 0 otherwise 

0=KBLK-0 c +(1-KBLK)(^ 6 


• Input from the main code 

1. Geometric quantities: z, y, X( , Xrh y ( , y v , J 

2. Flow variables: y, p, u, V, W, C w , C s 

3. Vectorization parameters 

4. Boundary parameters 

• Output 

1. To the main code: \i t , T { j 

2. For post-processing: K, e, y+, y n , fll 
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OVERVIEW OF PROBABILITY DENSITY FUNCTION (PDF) MODELING AT LeRC 

D.R. Reddy 

Internal Fluid Mechanics Division 
NASA Lewis Research Center 
Cleveland, Ohio 


OBJECTIVE 

Accurately model the effect of turbulence on 
chemical reactions in a fluid flow 


APPROACH 

Use Probability Density Function (PDF) model - 
Express dependent variables as functions 
representing statistically realizable events 


POSSIBLE MODELING STRATEGIES 

1 . Evolution PDF - solve for function 

a. Joint PDF for velocities and chemical species 

b. Joint PDF for only chemical species 
& energy 

2. Assumed PDF - function prescribed 
Limited range of applicability - 

reaction time < < or > > turbulence time scale 
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CURRENT APPROACH 


• Develop evolution PDF model for compressible 
reacting flows & extend to spray combustion 

• Solve for joint PDF for species and energy 
using Monte-Carlo technic|ue 

• Couple with conventional CFD codes 


AREAS OF IMPACT 

NOx Prediction - HSCT and AST application 


Spray combustion - swirling turb. reacting flows 


Scramjet flow path analysis 


Ignition kinetics - prediction of blow-off, etc. 


Combustion instability studies 


CODE FEATURES 


• Modular - can be coupled with any CFD code 

• Applicable for compressible flows with 
discontinuities 

• Monte-Carlo solver for generalized curvilinear 
coordinate system 

• Easily adaptable for parallel computation 
(currently under progress) 


CURRENT STATUS 

• 2-D and axisymmetric version released 

(default H2-air chemistry - 5 species) 

- parallel version to be released 

• 3-D version demonstrated for supersonic 
combustion (jet in cross flow) 

- validation planned for HSCT-type configurations 

• General chemistry (CHEMKIN) 

- Hydrocarbon spray combustion case currently 
under study 

• CFD codes used - RPLUS, ALLSPD, & SIMPLE-type 
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FUTURE PLANS 


• Further application/validation of 3-D model 

• Improved closure models - mixing and turbulence 
(use available DNS data) 

• Parallel processing - workstation clusters 

• Unsteady applications - long-term 

• Extend scope of impact 
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PDF METHODS FOR TURBULENT REACTIVE FLOWS N95. 27900 

Andrew T. Hsu 
NYMA, Inc. 

NASA Lewis Research Center 
Brook Park, Ohio 


OUTLINE 


• Motivation 

• PDF modeling of reactive flows 

• The Lewis PDF module 

• Validations and applications 

• Current research 

■ Technology transfer 


COMPUTATION OF TURBULENT COMBUSTION 


Conventional 
Model: 
Egregious 
• Results 



PDF Model: 

.••.ipThe cn uu 

Hope 





governing equations 


p,< + (pit;),i = 0 

+ (pUjU{)j = ~p,i + r i;j 
(p£).< + (pUjE)j = — <Ji,i + 'I’ 

(p%),i + (p«i^fc)j .= (p-O^j)j + 


A,, = 


dA 

dt 


Aj 


dA 

dxj 


CLOSURE PROBLEM; 


u t = u, + u \ , 

y t = K + y(, 


u"u'} 


YM 


Turbulence Modeling 

Analogy of shear stress: diffusion model. 


pi Wi 


??? 


But in general: 


pwi = pu>(Y\, Yli, ) 
•pWi # pi^Vi, ) 
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PDF Modeling of Turbulent Reactive Flows 
Current status 

. Assumed PDF (Spalding, 1971; Gosman & Lockwook, 

1973;...) 

0 Advantage: simple, fast. 

0 Disadvantages: Need unique mixture fraction; assumed 
shape may not be real. 

♦ Composition PDF (Pope, 1976; Dopazo &: O’B rian, 

1974) 

0 Advantage: Reaction rate created exactly; existing mo- 
ment closure codes easily adapted. 

o Disadvantages: Turbulent diffusion needs model. 

. Velocity-Composition joint PDF (Pope &c Chen 19S0. 

Pope 1981) 

0 Advantage: Reaction rate treated exactly; no dffusion 
model needed. 

0 Disadvantages: Models for velocity field relatively un- 
tried; Require more computer resource. 


PDF Modeling of Turbulent Reactive Flows 

• Objective: 

0 Develop models that can accurately simulate finite rate 
chemical reactions in turbulent flows. 

0 Develop and validate independent PDF modules. 

0 Technology transfer. 

• Criteria 

o Accuracy and robustness. 

0 Practical in terms of today's computing power. 

O Easy integration with existing industry computational 
platform. 
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PDF Modeling of Turbulent Reactive Flows 
. Approach: 

C> Joint pdf method for scalar compositions. 

0 Moment closure schemes for velocity field. 

0 Develop hybrid solver consisting of Monte Carlo method 
and finite-difference/finite- volume method. 


PDF Modeling of Turbulent Reactive Flows 
♦ Current status (Lewis) 

(pP),t -1- (p < Uj\Yi , h > P) j +{pwjP),y i 
= (. D t P-j)j + M\P) - (. S„P ). h ■ 


o Continuous mixing model developed. 

O Model for compressibility effect proposed. 

o 2D and 3D Monte Carlo PDF module developed. 


o 

0 


Validation studies. 

Code released l.o industry during a workshop. 
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Validation Cases 

Scalar field in homogeneous turbulence. 
Oblique shock. 

2D supersonic hydrogen combustor. 
Axisymmetric supersonic combustor. 
Piloted flame near extinction. 












Scalar field in homogenous turbulence 
pdf compared with Gaussian distribution 


Current model 



: 5 urc 2. A*r«pto*ie pdf dilution for * .uiir 
turbulence. — p<«« 1 r - 0<Jc, l ’ ’ * 


Modified curl's model 



- < -j > ;/u 


• U, u rv ! . /Va«tn|>'.oh: pdf distribution for « »c»Ur in homo 
turbulence. — tnodil'cd Curt ^«dc!; - - C*u*»'*J' 


Scalar field in homogenous turbulence 
3rd and 4th moments compared with 
Gaussian 


Current model Modified curl's model 



Evoluliou of nvoa*nls Lt;m the proem oodel. — 
iliidird delation. • • ■ 0.1 x fourth teslfjJ jrocseur. - - 
0.01 x tilth tciitril moment, ° 0.1 x fourth moment for 
distribution, O 3.01 x sixth moment for CUuuian 
dlj tributioo. 



7'irn? 


Involution of motwilj from the modified Curl 
model. - staoditJ dsvUlioa. ■ • • 0 01 x -otirtb «mr*! 

... 3c:c ^, -- 0.0001 x jixlh ccstfij honest, o 0. 01 x fourth 
wu.i.e3l'f« Ciu.fi^ distribution, C 0 -000 1 x filth moment 
for C*uui4Ji distribution. 
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Temperature across an oblique shock: 
pdf solution compared with analytical 
prediction. 


1.30 

1.25 

' □ pdf solution 

1.20 

— onolytico 

1 

1.15 

- 

O 

\ 1.10 
»— 

- 

0 

3 

1 .05 
i on 

c 

0 


i .UU 

0.95 



n on ■ 

' • t . 
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X 






Supersonic hydrogen combustor 
Mole fraction: 

pdf solution compared with exp. data 



Coaxial burner: geometry and test condition 


(Exp. Cheng, 

, et al. 

1991) 



([ Cnt Coodinctu 

0 ydro^cn 

jet 1 Outer Jet | 
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ii 
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Mean H20 mole fraction 
Coaxial burner 

pdf solution compared with exp. data 



r/0 * * * ' ' r / 0 


Piloted flame (Masri et al., 1994) 

Fuel: 45% CO, 15% H2, and 40% N2 
Flame close to extinction 



© 
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Piloted Flame 
Mean Temperature 







Current Projects 


0 Application of PDF module to emission predictions 
0 Incorporate general chemistry procedure. 

0 Incorporate spray models. 

0 Use parallel computing for the PDF module. 


Collaboration with industry and 
technology transfer 

. Features of independent pdf module: 

0 Easily coupled with any existing industry flow codes. 

0 Novel averaging scheme to reduce memory requiement. 
0 General chemistry package. 

0 Parallelized workstation version. 

• Technology transfer: workshops 

0 July, 1993; code released to 15 US institutions. 

0 October, 1994. 


257 



■I 


A COMPOSITION JOINT PDF METHOD FOR THE MODELING OF SPRAY FLAMES 

M.S. Raju N95- 27901 

Nyma, Inc. 

NASA Lewis Research Center 
Cleveland, Ohio 


INTRODUCTION 


• ™|COMPOSIT.ON JOINT PDF METHOD HAS 
BEEN USED TO MODEL A WIDE CLASS OF 
GASEOUS TURBULENT REACTIVE FLOWS 
(S.B. POPE) 

• NONLINEAR CHEMICAL REACTION RATES 
COULD BE EVALUATED WITHOUT ANY 
APPROXIMATION. 

• AN EXTENSION OF THE PDF METHOD TO THE 
MODELING OF SPRAY FLAMES. 

• EVALUATE THE LIMITATIONS AND 
CAPABILITIES OF THIS METHOD -IN THE 
MODELING OF GAS- TURBINE COMBUSTOR 


Composition Joint Pdf Transport Equation 


PP.t + PW.Pj, + [pw a (}0p]^ = 

{Mean convection} {Chemical reactions} 

-Ip < «? I ± > pU -Ip< - p J°^ I ± > 

{Turbulent convection } {Molecular mixing] 

~[p < I ± > pU 9 

{Liquid — phase exchange] 


• P 

• w a 


< t<" f 0 > 




1 ja 
•Sa | 0 > 


Density- weighted joint pdf. 
chemical source term for the a-th 
composition variable, 
conditional average of Favre velocity 
fluctuations. 

conditional average of scalar dissipation, 
conditional average of liquid-phase source 
term for the a-th composition variable. 
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Modeling Aspects of the Pdf Transport Equation 

• < u" [ ^ > is modeled using a gradient-diffusion model. 

• < I | t/j > is modeled using a variant of Curl s model. 

• The new term < ±s a \ tp > involving the conditional average 
of liquid-phase source term is modeled based on the average 
values of species and enthalpy: 

< [ ^ >= Y 

for <f> a = Ya>(X — 1^2, .. M s = <T — 1 

<~s a \lt> >= TT" Y n k m k{ — h.cII + hk, - 4> a) 

p — pt\v 

for <t> a = A. 


MODELING ASPECTS 


THE MODELED PDF TRANSPORT EQUATION 
PROVIDES THE SOLUTION FOR 
AND TEMPERATURE FIELDS WITH THE MEAN 
VELOCITY AND THE TURBULENT 

diffusivity and frequency PRoy‘9,1 R-Ai 

INPUTS FROM THE CFD SOLVER AND THE 
SPRAY SOURCE TERMS FROM THE 
LIQUID-PHASE SOLVER. 

THE MEAN FLOW AND TURBULENCE 
EQUATIONS ARE SOLVED BY A 
CONVENTIONAL CFD SOLVER WITH THE 
MEAN SPECIES AND TEMPERATURE FIELDS 
PROVIDED AS INPUTS FROM THE PDF 
SOLVER AND THE SPRAY SOURCE TERMS 
FROM THE LIQUID-PHASE SOLVER. 

THE LIQUID-PHASE EQUATIONS ARE 

FORMULATED IN .KStc 

COORDINATES WITH APPROPRIATE 
CONSIDERATION TAKEN INTO ACCOUNT OF 
THE EXCHANGES OF MASS, MOMENTUM, 
AND ENERGY BETWEEN THE TWO PHASES. 
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NUMERICAL METHOD 


• Mean-Flow and Turbulence Equations 

- Axisymmetric, Unsteady. 

- Incompressible Navier-Stokes (Variable-Density). 

- A Standard Two-Equation k-e Turbulence Model. 

- A Pressure-Based CFD Solver Based on the 
SIMPLE Algorithm of Patankar and Spalding. 

• Liquid-phase Equations 

- The Spray Model (Raju and Sirignano). 

- Dilute Spray Assumption. 

- The ODE’s for the Particle Size, Velocity, and 
Location are Solved Using a R-K Method. 

- The PDE’s for the Internal Droplet Distribution 
(Vortex Model) are Solved by an Implicit Method. 

- Droplet Regression Rate is Based on Either a 
Gas-Phase Boundary Layer-Analysis or 
Low-Reynolds Correlation. 


NUMERICAL METHOD 


• The PDF Transport Equation 

- A Fractional Step Monte-Carlo Method (Pope). 

- Spatial Transport, Molecular Mixing, Liquid-Phase 
Source Terms, and Chemical Kinetics are 
advanced in a Series of Sequential Steps. 

- Vectorization 


• Interaction Between the Two Phases 

- Interpolation of the Gas-Phase Properties at the 
Particle Location Using an Area-Weighted 
Averaging. 

- The Source Terms Evaluated at the Particle 
Location are redistributed among the surrounding 
Computational Nodes Using an Area-Weighted 
Averaging. 
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CHEMICAL KINETICS MODEL 


• IT IS BASED ON A SINGLE STEP GLOBAL 
MECHANISM OF WESTBROOK AND DRYER 
FOR N-DECANE/OXYGEN COMBUSTION. 


• THIS GLOBAL COMBUSTION MECHANISM 
WAS SHOWN TO PROVIDE ADEQUATE 
REPRESENTATION OF TEMPERATURE 
HISTORIES IN FLOWS NOT DOMINATED BY 
LONG IGNITION DELAY TIMES. 



2140 MM I 




a 

o 

o 




Geometry of the combustion chamber. 
(El Banhawy and Whitelaw) 
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EXPERIMENTAL DETAILS 


• The experimental data corresponds to 
the following inflow conditions: 

inflow temperature = 310 K, 
air mass flow rate = 355 kg/h, 
air/fuel ratio =20.17, 
swirl vane angle = 45 deg, 
swirl number =0.721. 

• The reported error in the 
measurements is about 10 to 15 % for 
the temperature and about 15% for the 
velocity. 


Details of Fuel Injection 


• A fuel nozzle of swirl-atomization type was used. 

• The liquid fuel injection is simulated by injecting 
a discretized parcel of liquid mass at the end of 
each Attn/ection 

• The droplet-size distribution is given by: 


— = 4.21 10 * 

n L D»\ D 32 



Droplet size Distribution 


• The initial droplet injection velocity corresponds 
to: u k = 11.0 m/s, w k = 6.1, and v k = 0.5 — 2.5. 
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PARAMETER SELECTIONS 


• The computations were performed on 
a grid with a mesh size of 60x60. 

• The PDF solution is obtained by 
making use of 250 particles per cell. 

• Dtg = Dtj n j ect j on = 1.5 ms, Dt^ = 

0.0375 ms, and DtMonte-Carlo = 
0.015 ms. 

• Two CPU seconds on a CRAY Y-MP 
per one Dtg and about 2 to 3 CPU 

hours - 4000 time steps - for the 
solution to reach steady state. 



Velocity vector plot 



Temperature contours and droplet locations. 


264 



Mean temperature, K Mean temperature, 



Axial d I s t dnc e . m 


Spray evaporation rate. 






Near wake radial profiles of temperature. 
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.52 







Near wake radial profiles of velocity. 



Schematic of an open spray burner. 
(Dan Bulzan of IFMD at LeRC) 
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CONCLUDING REMARKS 

• The comparisons show that the general features of 
the flowfield are correctly predicted by the present 
solution procedure. 

• The present solution appears to provide a better 
representation of the temperature field, particularly, in 
the reverse-velocity zone. 

• The overpredictions in the centerline velocity could be 
attributed to the following reasons: 

- The use of k-e turbulence model is known to be less 
precise in highly swirling flows. 

- The swirl number used here is reported to be 
estimated rather than measured. 
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IMPROVEMENTS AND NEW FEATURES IN THE PDF MODULE 


A.T. Norris N95- 27902 

Institute for Computational Mechanics in Propulsion 
NASA Lewis Research Center 
Cleveland, Ohio 


Overview 


• Modeling: What models are used in this package and what 
are their advantages and disadvantages. 

• Numerics: Describe how the PDF model is implemented 
and what are the features of the program. 

• Future Developments: What can be expected in the future 
from the NASA Lewis PDF code. 
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PDF Modeled Equations. 


Exact scalar PDF transport equation is: 






OXi — o<pa 


_ f) JCt 


dXi - 




(i> 


• Terms on the LHS exact - need to model the four terms 
on RHS, corresponding to turbulent convection, molecular 
mixing and the pressure term. 


Turbulent Convection 


• This term is modeled as a simple gradient diffusion process. 

(2) 

• D t is the turbulent diffusion coefficient, equal to the eddy 
viscosity. (Assume unity Schmidt) 

. Disadvantage: Counter-gradient diffusion known to occur 
in some pre-mixed flames. 
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Molecular Mixing. 


• Molecular mixing can be viewed as process which changes 
the shape of the scalar PDF without affecting the mean. 

• Molecular mixing is modeled by two models: A coalescence/ 
dispersion model (Hsu and Chen) and a relax-to-mean model 
(Dopazo). 

• Advantages of both models is that they are simple and 
readily adaptable to any number of scalars. 

• Disadvantages are the relative lack of physics in the models. 


Pressure Term. 


• Pressure term model is based on second order closure mod- 
els for compressible flows (eg. Sarkar). 




m 

dt 


+ m£-+ o.8„ w ^ 

-I- 0.15pP r Mt — 0.2peM? 


(3) 


• Advantages are that model is tried and tested in finite 
volume codes. Disadvantage is that only the mean pressure 
can be used for model. Idealy we would like a stochastic 
process for two state variables. 
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Numerics 


• Solution of scalar PDF transport equation achieved by a 
particle based Monte Carlo scheme. 

• PDF represented by an ensemble of particles, each with a 
composition and enthalpy. 

• PDF evolves by the motion of these particles in physical, 
scalar and enthalpy space, by exact and modeled processes, 
eg. Convection, reaction, mixing. 

• Statistics (eg. means) obtained by averaging over ensemble 
of particles. 


Numerical Details - Monte Carlo Scheme 


• Module based on cell-centered quantities. 

• PDF method is a nodal one. ie. All particles reside at the 
center of the cell, and can move only to neighboring ceils. 

• Number of particles at each node is the same and the 
number remains fixed. 


• Statistics obtained by averaging over particles at each node, 
and also by time-3veraging. 
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Numerics - Convection 


• Upwind scheme used for mean convection, and central dif- 
ference for turbulent convection. 

• Evolution achieved by moving particles in from adjacent 
nodes. Particles are selected at random. 

• Fractions of particles are treated by random convection. 


if 6.3 particles then 


6.0 70% of the time 

7.0 30% of the time 


( 4 ) 


Numerics - Reaction 

• Although reaction source term treated exactly, several dif- 
ferent numerical schemes are needed. Timing figures are for 
% of time spent in the PDF part of the code on SPARC II 
workstation. 

• No reaction: For scalar mixing calculations. 

(Timing: 17.2 %) 

• Equilibrium reaction: Assume reaction proceeds at infinite 
speed. Table of equilibrium composition as a function of 
mixture fraction obtained from seperate CHEMKIN routine. 
(Timing: 24.4 %) 
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Numerics - Reaction Cont. 


• One-step global reaction schemes. Westbrook and Dryer 
global reactions integrated for each time-step. 

(Timing: 51.1 %) 

• Tabulated reaction incriments. Users create their own table 
of composition incriments as a function of scalars using the 
adaptive tabulation scheme provided, plus the users favourite 
reduced mechanism. 

(Timing: 58.9 %) 

• Chemkin full mechanism integration. Very slow and not 
recomended except for parallel applications. 

(Timing: 97.8 %) 


Numerics - Averaging 

• To reduce statistical error in the evaluation of the mean 
scalar quantities (without increasing the number of particles 
per node), time averaging is employed. 

• A weighted time average is used to give more weight to 
recient values and less to those in the far past. 


(( <p))t = -■ x r (Wt + w t{(<i>))n- 1) 

w t + 1 

Wt = Ct(w t - 1 + 1 ) 


( 5 ) 

( 6 ) 
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Numerics: Misc. 


• A portable random number generator is now included in the 
module, set up for 32 bit machines. 

• A time step check is now performed to ensure boundedness 
of the PDF solution, ie. no negative numbers of particles. 

• Rplus/PDF release ported to workstation enviroment. K- 
epsilon now standard turbulence model. 


Future Work. 

• Release of 3D version with new improvements. 

• Implimentation of parallel processing for distributed cluster 
environment. (PVM based) 

• Include model for another state variable to close PDF mod- 
eling. 
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